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ABSTRACT 


A system is proposed for grid allocation and differencing of apparently general applicability to purely marching- 


type systems of equations of fluid dynamics. 


The method is based on casting of the equations into the conservation 


form, which then permits use of a staggered space-time grid system with interpolations required only in certain 


linear terms. 
experiments have been successfully performed. 


The method is illustrated by application to two systems of equations, on one of which numerical 
Advantages and drawbacks of the method are 


described in com- 


parison to other currently used grid systems, and the possibility and desirability of parametric simulation of turbulent 


eddy exchange processes are discussed. 


1. INTRODUCTION 


The recent development of methods of numerical inte- 
gration of the equations of fluid dynamics in a more-or-less 
primitive form, i.e., of first order, except for diffusion 
terms, is focusing attention upon new aspects of numerical 
analysis. During the early development of geostrophic 
or balanced models which involved a vorticity equation, 
perhaps the most pressing practical problems were to 
solve rapidly and accurately certain rather complex 
elliptic second order partial differential equations. In 
dealing with undifferentiated systems one finds that the 
basic equations are much simpler to apply but must, for 
the sake of computational stability, be applied at very 
short time intervals compared to the time scales of the 
significant meteorological phenomena. This is because 
the normal Courant-Friedrichs-Lewy stability criterion 
requires essentially that motions or waves allowed in the 
system not be able to travel from one grid point to the 
hext in one time step, and the undifferentiated systems 
usually allow faster-moving waves than those of meteoro- 
logical interest. It is also noted that the boundary con- 
ditions and various consistency requirements, if not more 


complex, are more critical when a single computation run 


may include several thousand time steps. Under these 
circumstances it is desirable to examine the time-space 
grid data allocation and differencing scheme with a view 
toward eliminating, if feasible, some of the redundant 
time resolution and thus save time and computation 
expense. It also appears that when integrating a set of 
non-linear equations over a moderate number of time 
steps, say several hundred, a certain form of instability 
arises, Which is related to spatial truncation error in the 
non-linear terms. This non-linear instability, discussed 
by Phillips [7], seems to occur much more rapidly in 
primitive equations models than in those using a vorticity 
equation, though it is uncertain whether this is due to 
differences in the physical or mathematical behavior of 
the systems. At 
possible devices for elimination of this instability, which 


any rate there appear to be several 
will be mentioned later. 

As a method of eliminating some of the computational 
redundancy, Eliassen [3] proposed au method of handling 
a two-level primitive equations model, in which variables 
are staggered in space and time in a somewhat complex 
manner. The principle on which the system was based 
was that the linear terms of the equations would be 
available at the correct grid points with minimal trunca- 
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certain non-linear terms would be 
formed from The 
reason for doing this, presumably, was that the linear 


that 
products of 


tion error, but 
interpolated values. 
terms of the motion equations, the pressure gradient and 
Coriolis terms, are generally an order of magnitude 
larger than the others and errors in these terms might be 
thought to be most damaging. On the other hand the 
advective terms may contribute as much, and in many 
cases much more, toward significant rates of change of 
quantities than the combination of the nearly balanced 
As long as truncation error does not disturb 
the linear terms 


which is most pronounced for the larger scales of motion, 


linear terms. 
the existing quasi-balance between 


there seems no reason to require interpolation to be done 
only on the advective terms. It will be shown that the 
truncation error of interpolation in the proposed system is, 
at most, of the same order as that arising from the finite 
difference pressure gradient term, but that the latter is 
evaluated with less truncation error in the Eliassen grid 
than in that proposed. 

Hinkelmann [4] and Smagorinsky [13] have successfully 
applied hydrostatically filtered equations of two-dimen- 
sional flow without the use of time or space staggering of 
the both of the external 
gravitational motions were filtered out by the vanishing 


variables. In these models 


boundary conditions on dp/dt, and it was therefore neces- 
Both 


used the marching equations in the conservation form, and 


sary to solve a Poisson equation at each time step. 


it was found that non-linear computational instability 
ensued, in the absence of specific damping terms, within 
two or integration. 
Phillips [8] has recently developed a barotropic divergent 


three days of commencement. of 


(free surface) primitive equations model, in which the 
Kliassen grid system has been applied, with central time 
and space differencing, and also with a one-sided space 
difference scheme |9] oriented according to the wind 
direction and related to that originally proposed by 

The latter 
damping 


Courant, Isaacson, and Rees [1]. method 


introduces some computational and seems 


effectively to prevent development of non-linear 


instability. 


2. PRINCIPLES OF PROPOSED SYSTEM AND 
APPLICATION 


A method will now be presented for stepwise integration 
of initial-value boundary-value problems using non-linear 
hydrodynamic equations. The 
generally applicable to meteorological problems when the 


method appears to be 
system of equations is of an explicit marching type, that 
is When there are no physical approximations or con- 
straints which involve the solution of an elliptic equation 
ateach time step. It is evidently not particularly suitable 
for application to, for example, the barotropic vorticity 
equation, or to primitive equations models such as those 
of Smagorinsky and Hinkelmann in which’ external 
gravity Waves are filtered out, because in these systems 


A principal 


further interpolations would be necessary. 
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feature of the method is the casting of the equations into 
the conservation form, the form in which all advective 
terms appear as flux divergences. When this is done it 
becomes very natural to set up a central time and space 
differencing system on a space-time staggered grid, in 
which interpolations may be required only in certain 
linear terms. The method of interpolation is largely 
optional, but time interpolation is generally most accurate, 
If frictional terms are present (and it is believed that they 
generally should be) these are most naturally computed by 
forward differencing, which then simultaneously preserves 
linear computational stability in terms. The 
method will be illustrated by application to a free surface 
model similar to that of Phillips. 

If the atmosphere is considered as a homogeneous in- 


these 


compressible fluid with a free upper surface, whose height 
is @/g, the hvydrostatically _ filtered 
equations of motion may be written, in Cartesian coordi- 


two-dimensional, 


nates, as 


na ou ~ nee 1 Or 
Or or ‘Or, o OF 
where €;, is the permutation tensor, equal to plus or 
minus unity for 7, 7, k=1, 2, 3 or 2, 1, 3 respectively. 
Upon application of the 
without assuming constant viscosity, the frictional stress 


viscosity hypothesis, but 


may be written as 


- (Ou; , OU; 26,,; OU, 

ok ( + —_—— ): 9 
Or, OF; 44; OF, 

The repeated index implies summation, so that 6,, equals 

two in this 

system whose characteristic Reynolds number was small 


case. If we were considering a_ physical 
for example a “‘dishpan”’ experiment, we could reasonably 
assume @K constant and the stress terms would reduce te 
For the atmospheric 


system Wwe want the stress tensor to describe the transfers 


the common Laplacian form. 
of energy, momentum, ete., performed by turbulent eddy 
motions, the scales of which range between that of the 
smallest explicitly described motions and that of molecular 
dissipation, [14] proposes that for fully 


turbulent flow A be made proportional to the deforma- 


Smagorinsky 


tion, 1.e. 


K (kl) OK 5) 


where / is the seale of the smallest resolvable motions (the 
grid scale in a finite difference formulation) and f/ is a 
universal constant of order unity. From consideration ol 
the similarity of the above expression to that frequently 
applied in turbulent boundary layer theory we consider é 
to be analogous to the Karman constant. The appear- 
ance of a particular length seale in the definition ts a 
reflection of the essential correspondence of the frictional 
stress to the Reynolds stress, also defined in terms of the 
same length scale. 

The form proposed by Smagorinsky is similar to one 
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applied by von Neumann and Richtmyer [10] in their 
“pseudo-viscosity”’ method for representing shock waves. 
It has been introduced here because of its seeming impor- 
tance in elimination of the non-linear computational 
instability which is otherwise permitted by the proposed 
grid differencing scheme (and many others in current use). 
The physical significance and justification of these terms 
will be shown by Smagorinsky and are not discussed 
further here. 
The continuity equation is a prediction equation for the 
height of the free surface, written as: 
Op O( du ;) 


a 


=()- ( 
ot Or; 4) 


If (4) is multiplied by u; and combined with (1), a momen- 
tum equation is obtained, that is 
O(¢?/2) Ori; 


(pu ,u,) — fe; 36Uj+ . (5) 


O(du;) , O 
- or Or, 


Ot 


We now assume that motions within a rectangular area 
of dimensions L, are cyclicaily symmetric in all directions, 
so that the boundary conditions may be written, in terms 
of an arbitrary dependent variable x, as 


x (41, t2t Ly, t). 6 


Equations (4), (5), and (6) are a complete set of prediction 
equations for the dependent variables, @u;, @, written in 
the conservation form. 

We divide the 1, x L, area into pq grid squares, each of 
Letting / 
and m represent the index numbers of grid intersections, 


area A?, where p and qg are both even numbers. 


and n a time index, we define 


WP — XA: MA, NAHE, HI (7) 
the latter two equalities proceeding directly from (6). 
Let us now specify all the dependent variables at the 
even intersections (/+m even) for even n and also all of 
We then may 
by finite difference equations as 


them at the odd intersections for odd n. 


approximate (4) and (5 


follows: 


Af 
@ oO - (Ou, ) (pu) . 
(OUs) 1 m4 (Puls), » ‘ Ss) 
r ! r Af : 9 
ou (ou,):. — A (OU)° + O°/£— 711) 141.7 
(@U,*-+ O°/2— 711) »-1. mt (OU) Ue2—T 12) 1. m3 
(PU U2—T12) 1, m—-1) " +2fAt(g@us)(", (9) 
— ae 
oO (OUs); | (QU; Us 791) 14+1.m 
A 
(PU U2— 721) 1-1, mt (PU2? + G?/2— To2) 1, m+ — (PUe* 


+? /2— 92) 1, m—1)"" — 2fAt(ou)(”%. (10) 


MONTHLY WEATHER REVIEW 61 


The stresses and the viscosity coefficient A may also be 
approximated, from (2) and (3) as: 


(ri) f= — (722) (= (0), Ai, /A)[ (tr) p41, m 
(21) ¢—1, ms (Ue) 2. m+it (Ue) 2. m—1)" (11) 
(Tio) 0 m= (Tardem (o)"),A me /A)[ (tr) &, tt 
— (Uy) 2. m—1-t+ (Ue) 141. m— (Ue) 1-1. |" (12) 
Ky m=FPA{ [ (Ur) p41. m— (U1) 2-1, m— (Ue) 2, mg + (Ue) 1 m—1 |? 


+ [ (1) 2, m+1— (U1) 2, m (2) 141, m— (U2) 1-1. m)*} 4/7. (13) 
The bar over the Coriolis terms in (9) and (10) indicates 
that some form of interpolation or replacement must be 
done, since the quantities involved are not defined at the 
proper time-space grid points. 

Four possible methods of interpolation are suggested. 
Spatial interpolation would involve replacement of the 
barred quantities in (9) and (10) as follows, according to 
either a 2-point or 4-point formula. 


(pu;)/") is replaced by 1/2 [(@ti)), m41 + (OU) 1, m—1] 





and (14) 
(@uz)/", is replaced by 1/2 |(@u.) (Oto) 7-1. m] 
or 
(pu;)/". is replaced by 1/4 [(@u,) 41, 
(bu; );—1 mt (OU;) (pu,), , . (15) 


Other methods of spatial interpolation would involve 
somewhat larger truncation errors. For time interpola- 
tion one would use similar replacements us follows: 
(ou,)/". is replaced by 1/2 |(@u,);" (pu,)\"—" (16) 
where (9) and (10) would then have to be linearly com- 


bined to eliminate the implicit terms. The resulting 
system will be shown to have negligible truncation error 
in the Coriolis terms, as well as a slightly less restrictive 
computational stability criterion. A simpler method in- 
volving somewhat larger but still acceptable truncation 
errors is a forward-backward scheme similar to that in- 
troduced by Courant, Friedrichs, and Lewy |2| for a linear 
wave problem. In this case the barred quantities in (9) 


and (10) are replaced as follows: 


(@u;)}",, is replaced by (@u,)/"," >) 
(17) 
(@uy)/") Is replaced by (@us);", 
It may be noted that the frictional terms in (11), (12 


and (13) involve velocity derivatives at time index (n—1) 
weighted by @ values at time index (nr). This inconsist- 


ency (if such it is) seems from experience to be totally 





62 
unimportant. It has, in fact, been convenient to compute 
similar frictional terms over considerably longer time steps 
than the remainder of the calculations, and little error is 
so introduced. It is well known, however, that straight- 
forward central differencing for the diffusion equations 
may lead to linear computational instability (see, e.g., 
Richtmyer [10}). 

With regard to initial conditions, it may be desired to 
begin computations with values of variables all pertaining 
to the same time, rather than in the staggered grid. In 
this case it is satisfactory to make an initial half-time step, 
and similarly to interpolate between time steps in order 
to obtain a simultaneous display at a later time. 


3. TRUNCATION ERROR AND COMPUTATIONAL 
STABILITY 


We shall not attempt to give a complete discussion of 
the computational stability and truncation error of finite 
difference formulations of (1)-(4). We have, however, 
determined the stability properties of linearized versions 
of several finite difference equivalents of the frictionless 
system and first approximations to the truncation errors 
of various terms in the motion equations. This material 
is assembled in tables 1 and 2. The expressions given for 
truncation error are the lowest order terms of the Taylor 
series, and are presented in this form for comparative 
purposes rather than for any attempt at numerical eval- 
uation. The grid formulations used include (a) the nor- 
mal advective form, based on equation (1) with central 
differencing and no interpolations; (b) the proposed 
staggered grid system with 2-point and 4-point space in- 
terpolations, time interpolation, and forward-backward 
treatment of the Coriolis terms; and (c) the Eliassen grid 
as used by Phillips [8]. Figures 1 and 2 illustrate the 
essential difference between these grid systems. In the 
normal advective scheme all variables are defined at all 
intersections of figure 1, while in the proposed staggered 
system all quantities are defined for even time-steps at 
grid points marked “FE” and for odd time at points “O”’. 
In the Eliassen scheme the geopotential, @, is defined on 
grid ‘“‘A”’ of figure 2 for even time and grid “D” for odd 
time, uv, is defined on grid “B” for odd time and grid “C” 
for even, and w on grid “C”’ for odd time and grid “B”’ 
for even. The grid interval shown on figure 2 was altered 
to (1/y2) of that defined by Phillips in order that the 
truncation errors for the various systems be comparable. 
Thus in all cases the total number of points describing a 
given field, combining odd and even time steps, is given 
by the ratio of the total area to A’. 

Table 1 shows that the non-viscous linear computational 
stability characteristics are similar for all the grid-differ- 
encing systems considered. The differences lie in’ the 
Coriolis term, which is ordinarily two to three orders of 
magnitude smaller than that arising from the external 


gravity wave propagation. It may be easily demon- 
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Figure 1.—Grid system for normal advective scheme, in which all 
variables are defined at all intersections; and for proposed 
staggered system, in which all quantities are defined for even 
time-steps at points E and for odd time at points O. 


strated that the friction terms, computed in the form de- 
scribed, exert an effect on the stability criterion of at 
most the same order of magnitude as the advection. 
Upon careful examination of table 2 we find that 
truncation errors also do not present any clear superiority 
of one system over another. Truncation errors in the 
geostrophic terms are certainly minimized in the Eliassen 
scheme, because of the smaller differencing interval in the 
pressure term. Nevertheless the unavoidable error of 
first-order differencing of the pressure term is of the same 
order as that of spatial interpolation of the balancing 
Coriolis term, so that introduction of the latter breeds no 


Table 1. Linear com putational stability criteria 
" , A 
Normal advective TF \2o+(fA 
Al 
. . A ; 3 
Space-interpolated staggered Ar uj +y2o+ (fA)? 
ee A 
rime-interpolated staggered uj +y2o 
Al 
: A ; 
Forward-backward staggered Al ui; +v2e@+ (fA)? 
sia A : 
Eliassen u;|+yv2o+ (fA)*2;2 
Al 
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FicurE 2.—-The Eliassen scheme in which the geopotential @ is 
A for even grid D for odd time, 

is defined on grid B for odd time and on grid C for even time, 
ind 42 


ne 


lefined on grid time and on 


is defined on grid C for odd time and on grid B for even 


new sources of error. Proper comparison of the various 
methods of interpolating the Coriolis term is dependent 
upon specification of characteristic time and space scales 
of motion. For typical meteorological motions (small 
gravitv-wave amplitude) the time interpolation has the 
smallest error, only a little larger than that of the time 
derivative The the 


scheme is considerably larger and may be nearly as large 


term. error of forward-backward 
as that of the space-interpolation method, depending on 
the spatial smoothness of the fields. Comparison of the 
truncation errors associated with the non-linear terms is 
more difficult. It that terms of all 


systems are of roughly the same size except the Eliassen 


appears, however, 
component along the wind, which is about four times 
larger 

All results available to the author indicate that all the 
systems considered above exhibit non-linear computational 
instability unless their short-wave components are effec- 
tively damped. Our understanding of this phenomenon 
is, however, severely restricted by the present lack of 
general methods of its analysis. Phillips’ [7] method may 
be used to prove instability, in some cases, but it can 
hever prove stability. Trial-and-error methods are rather 
ineflicient, with the large number of possible methods to 
choose among. Shuman [11] has, however, recently per- 
lorined one-dimensional numerical experiments, the re- 


sult. of which raise the possibi’ity that certain reasonable 
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TABLE 2.—Truncation errors of terms in equations of motion 


. . Op 
1. Pressure derivative term; e.g. 
Or; 
. P A’? 0%@ 
Normal advective ~ 
6 Ox,3 
St red oe 
‘ a rere 
8s ce) 12 Or,” 
— A? O°'o 
Eliassen . 
12 Or,' 


2. Coriolis term; e.g. fu. 


None 


Normal advective 


‘ , : f A? 0? (ou 
2-point space-interpolated staggered 


os Or; 
f A? 

{-point space-interpolated staggered ig © (ou 
ta) ) 


Time-interpolated staggered 


QD 2 Of? 
. / Ol ou 
Forward-backward staggered Al 
ta) Ol 
Eliassen None 
2 Th . Ou 
3. Time derivative, e.g. 
ol 
’ . At)? Ou 
Normal advective = 
6 ol’ 
‘ 1 ( At)? 08 ( pu) 
Staggered £ 
o 6 ols 
— Al Oru 
Eliassen 
6 Ul 


Ji 
4. Advection term along wind component; e.g. u 
Ou 
° » Oru 
Normal advective ‘ 
b Ul 
‘ l A’ O° Qu 
Staggered - 
o 6 Ol 
: A 2) 0 Ou 
Eliassen ( ; ) 
6 \ Ou Ou Ou 
_ Ou 
5. Advective term across wind component; e.g. u 
Ou 


. : a Oru 
Normal advective h.. 
b Ul 
. 1 A? OF ( @uyu 
Staggered 
o 0 oF 


Eliassen 
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difference schemes may be stable without damping. — If 
these results can be shown to apply to a more general 
model we may have a close counterpart to the spectral 
method which in some cases may conserve energy identi- 
cally. 

The eddy-viscosity term proposed by Smagorinsky 
obviously exerts a damping which is highly selective in 
scale. For a given amplitude of disturbance kinetic 
energy and a given grid interval the value of A is propor- 
tional to the finite difference approximation of the de- 
formation amplitude. When central differencing is used 
each contribution to this amplitude will be proportional 
to sin (27A/d), where \ is the wavelength of a motion 
component. Thus damping is maximized for wavelengths 
at and near 44 and actually vanishes at A=-2A. 
components with wavelengths between 24 and 44 are 


Since 


most directly involved in non-linear instability it is there- 
fore suggested that the viscosity term prevents the insta- 
bility indirectly. That is, it reduces the amplitudes of 
the intermediate scale (stable components which would 
otherwise generate small-scale (unstable) components. 
If this is the case we may have reasonable confidence that 
the important effects of the analytic viscosity terms (2), 
(3) are fairly well approximated by their finite difference 
formulation, 

Other methods that have been used to prevent non- 
linear instability do more or less violence to the pliysies. 
Probably the worst method is to use a constant viscosity, 
or its equivalent in a difference system or smoothing 
filter, 
effective the viscosity must be large enough so that the 


Experiments have shown that in order to be 


greatest characteristic grid Reynolds number of the field is 


of order unity, where by contrast we note that 
Smagorinsky'’s eddy viscosity is of such a form that the 
grid Reynolds number is everywhere of order unity. In 
order to simulate motions of a turbulent fluid one would 
then have to carry millions of grid points, otherwise the 
system would be essentially laminar and the large-scale 
components severel\ oversmoothed, Lax [6] suggested a 
space-time staggered grid for integration of equations in 
the conservation form. Forward-time differencing was 
used and variables at the previous time were space- 
averaged before differencing. This process effectively 
introduces a very large constant computational viscosity, 
essentially equal to A?/4 Af, and the grid Reynolds number 
is of order u/yo<1. On somewhat the opposite extreme 
of selectivity is the method used by Phillips [7] in’ his 
general circulation calculations, which consisted of Fourier- 
analyzing the motion field and removing all components 
with wavelengths less than 44. Somewhere between these 
extremes lies the one-sided difference scheme first: sug- 
gested by Courant, Isaacson, and Rees [1], in which the 
advective form of the veetor equation of motion is used 
and spatial derivatives are taken along the characteristics 
(streamlines) upward of the central point. Variants of 
this method have been used by Kasahara 15] and Phillips 
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[9]. Motions are damped, but rather selectively, and the 
general results appear to be very similar to those obtained 
using Smagorinsky’s eddy friction form, although the 
physical significance is unclear. 


4. OTHER COMPUTATIONAL FEATURES 


An advantage of the proposed system, in comparison 
with, for example, the Eliassen grid, is the vanishing trun- 
cation error of volume integrals (numerical sums) of the 
prediction variables. Smagorinsky [12] has shown that 
this is assured by application of certain auxiliary boundary 
conditions (mainly symmetry conditions) and these in 
turn prevent development of another type of computa- 
tional instability. The importance of this feature probably 
dependssomewhat on the use to which the particular predic- 
tion system is put. For general circulation studies, where 
integrated momentum budgets, heat budgets, ete. are 
of fundamental interest, it seems very desirable to know 
that, for example, the mean geopotential is only a function 
of the boundary conditions. For operational forecasting 
this may be of less direct interest, though perhaps com- 
forting knowledge. In any case it has proven to be 
valuable for “debugging” machine program logic, for deter- 
mination of machine errors, and for evaluation of round- 
off errors 

The proposed system has certain other practical advan- 
tuges for application to a high-speed electronic computer, 
not all of which are shared by other systems. First, only 
one set of variables need be stored for each point, in 
distinction to the two or three required for a non-staggered 
grid using forward-differenced viscous terms. Second, 
there are just two kinds of points, rather than four as in 
the Eliassen system, which makes the logic somewhat 
simpler. On the other hand, application of any staggered 
grid system adds definite logical complications to a 
machine program, thus increasing the programming and 
check-out time and, to a small extent, the time spent by 
the machine in logical testing. 


5. APPLICATION TO OTHER MODELS 


The proposed method was originally developed for, and 
applied to, integration of a set of two-dimensional (2, ts) 
equations used for the simulation of dry convective 


motions. The system used was the following: 


O(pu,;) re) Op 
+ ( 
or TS eu y)" or + 908:3 or 
e+ (pu,)=0 
I 
O( pA) ra) ofl 
(pu 6) 
Ot or or 


where pressure is obtained from the equation of state in 
the form 
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R pd " 
7 


pP 
m4 
and the frietion-diffusion terms were of a similar form, 
though more complex, than those described above. 

The grid scheme outlined above is clearly applicable to 
this set, where the linear density term in the vertical 
momentum equation must be interpolated. In this case 
it is possible to interpolate density either spatially or time- 
wist Integrations have been successfully performed 
using two-point (vertical) and four-point spatial interpo- 
lation and linear time interpolation, the results of which 
will be reported elsewhere. 

If the above set of equations were modified to exclude 
sound wave propagation, by the time 


derivative of density and formation of a vorticity equa- 


elimination of 
tion, it would then not be practical to apply the proposed 
system. On the other hand, it can evidently be applied 
to a multi-level hydrostatically filtered baroclinic model, 
provided that the vertically integrated divergence is not 
constrained. The applicability of the system depends not 
on the scale or complexity of the model but only on the 
absence of a constraining elliptic differential equation. 


6. SUMMARY 


It has been demonstrated theoretically and, in one 
particular case, practically, that the staggered grid system 
here proposed can be advantageously applied in the inte- 
gration of purely marching type hydrodynamic-thermo- 
dynamic models. The method is of a similar type to 
Eliassen’s but essentially opposite in its basic principles, 
in that here the linear terms may be interpolated, spatially 
or perhaps preferably in time, while the advective terms 
are applied with centered differencing in the conservation 
form 
energy from grid-scale motions by parametric simulation of 


Physical considerations suggest the abstraction of 


eddy exchange processes, and Smagorinsky’s suggested 
method for doing this seems to be sufficient for mainte- 
nance of computational stabilitv. The essential features 
of the system have been described within the framework 
of a rotating incompressible barotropic fluid model with a 
free surface, while its practical application has been to a 
somewhat more complex model of a compressible fluid 


with thermal convection. 


to 
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DIAGNOSIS OF DIVERGENCE IN A THREE-PARAMETER 
NUMERICAL PREDICTION MODEL 


A. WIIN-NIELSEN 


ABSTRACT 


The equation for the vertical velocity in a quasi-non-divergent, three-parameter model has been solved for a 


certain simple flow pattern 


discussed from the solutions 


computed and compared with those obtained by other investigators 


The importance of vertical variation of static stability and of the horizontal wind is 


Examples showing the distribution of divergence relative to the synoptic systems are 


A general discussion of the factors influencing 


the mid-tropospheric divergence follows in section 4, and section 5 contains finally some remarks on the divergence 


in very long waves 


1. INTRODUCTION 


The distribution of convergence and divergence and 
r re 


of vertical velocity relative to the atmospheric flow 
pattern has been a very important and most intriguing 
problem in synoptic and dynamic meteorology for many 
vears. Due to the fact that these quantities have to be 


obtained by indirect methods numerous rules relating 
the distribution of divergence to the synoptic flow pat- 
tern have been formulated. Intensive synoptic studies 
Fleagle, [2]) have resulted in characteristic distributions 


Theo- 


retical studies by Charney [1] and others using a very 


of divergence relative to the troughs and ridges 


general continuous model gave as a_ by-product the 


divergence fields in baroclinic waves. General agreement 
between the results of the svnoptic and theoretical studies 
Was apparent. 

The different dynamical models applied in numerical 
weather prediction and in studies of the general circula- 
tion, except the non-divergent model, contain implicitly 
distributions of divergence and vertical velocity, but 
because most of the time integrations have been made 
vith one- or two-parameter models it has not been too 
teresting to study the predicted distributions since the 
vertical variations are constrained to very simple patterns 
like parabolic distributions of vertical velocity with 
pressure and linear variations of divergence 

The experiments with quasi-non-divergent two-param- 
ter models have shown that these models are not 
sufficiently 


strong baroclinic nature. 


accurate to predict the developments of 
Several reasons may be men- 
tioned for these somewhat discouraging results, but it is 
evident that the two-parameter model which assumes a 
constant direction of the thermal wind is not able to 
describe the vertical variation of the temperature advec- 
lion pattern. Quite frequently we find regions of cold 
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air advection superimposed on warm air advection in the 
same vertical column and vice versa. The only way to 
incorporate such features is to introduce more information 
levels to get a better verticat resolution of the model 
From the operational point of view it is obvious that 
as few levels as possible are desirable. Not only does 
this decrease the time required to compute the forecasts, 
but the analysis problem becomes also less time consuming 
One may therefore ask whether the addition of one extra 
information level will change the behavior of the model 
to such an extent that we can expect substantial differ- 
ences between the three- and the two-parameter models, 
It is the purpose of this study to show that the distri- 
butions of divergence and_ vertical velocity which are 


inherent in the three-parameter model show a great 
similarity to the distributions obtained by Charney [1] 
Fleagle [2] and more recently Hinkelmann [3]. A treat 
ment of the three-parameter model as compared to the 
continuous case has the advantage that the solutions for 
divergence and vertical velocity considering relatively 
simple flow patterns can be expressed in such a form that 
the importance of the different factors contained in the 


We shall thus be 


able to consider the lportance of the vertical variation 


solution can be easily investigated 
of stability and the variation of wind speed with height, 
especially a crude measure of the curvature of the wind 
profile. Another advantage in this case is that we shall 
be able to study the solution for different horizontal scales 

The vertical variation of static stability will in this 
study be prescribed as a function of pressure in such a 
wav that we obtain the major part of its systemats 
climatological variation. It should be stressed that this 
variation does not necessarily correspond to the svnopti 
variability. 

It is necessary to restrict the investigation in this papel 


to flow patterns where the zonal wind is constant in each 
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isobaric surface, but of course can vary with pressure, 
and to waves of the simple sinusoidal tvpe. The results 
which can be obtained in this way are of course not directly 
applicable to the real atmosphere with its much more 
complicated flow pattern. Solutions for observed atmos- 
numerical 


pheric flow patterns can be obtained by 


methods. The diagnosis of the real atmosphere using 
equations very similar to those in this paper is presented 


by Cressmian in another paper in this issue. 


THE COMPUTATION OF VERTICAL VELOCITY 
IN THE THREE-PARAMETER CASE 


The three-parameter representation of the atmosphere 
allows a solution for the vertical velocity at two internal 
levels in each vertical column in addition to the two values 
obtained from the boundary conditions at the top of the 
atmosphere and at the ground. Using a finite difference 
form of the continuity equation we can consequently 
obtain three values of the divergence at the intermediate 
levels In this paper we shall for simplicity divide the 
atmosphere ito SIN lavers, each corresponding to p 6 eb., 
where Py, is the pressure at the ground (fig. 1 We shall 
assume that p 100 eb. and thus neglect the effect of 

The boundary conditions for the vertical 
dp dt, will be w 


is quite likely that it will be an advantage to assume w= 0 


topography 
velocity, a 0 for p=0 and p=p It 
at some level in the stratosphere in the practical applica- 
tion and also to include the effects of the topography and 
friction at the lower boundary 

The equation for the vertical velocity is obtained from 


the vorticity equation in the form: § 


Oc . Ow 
tV.9(6+-f=f 2.1 
ot Op 


and the adiabatic equation 


Ww 0. 6.4 


0 (OY oy o 
st\ap)t V-v( Op a 


In (2.1) and (2.2) W=k XVy is the horizontal wind 
assumed to be nondivergent, k au vertical unit vector, y the 
streamfunetion, ¢=Vey the vertical component of the 
dp/dt the 
adln @/Op a measure of static 


relative vorticity, / the Coriolis parameter, w 
vertical velocity, and @ 
stability ‘ais specific volume and @ potential temperature. 
The vertical derivative of the geopotential has been 
replaced by the vertical derivative of the streamfunction 
in the adiabatic equation (2.2) using the approximate 
relation (Phillips [5 


OY Op=1/f,.0¢/Op 2.0 


where @=gz is the geopotential, g the acceleration of 
gravity, and 2 the height of the isobaric surface. 


From (2.1) and (2.2) we obtain the w-equation 


9 e O*w oO ° ow 
oVw-+-f,* = Tn (V.0T (C+) “ e « (2.4 
Si Sp3 [ overt ) vi(ver oe) | ) 
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In the derivation of (2.4) it has further been assumed 
that the horizontal variation of the static stability can 
be neglected, while a vertical variation is still possible. 

The procedure to be followed is now to apply (2.4) at 
the levels 2 and 4 (see fig. 1), to approximate vertical 
derivatives by finite differences in a straightforward Way, 
and in this way obtain two equations both containing w, 
and w, as unknowns. The streamfunction, y, has then 
to be known at the levels 1,3, and 5. It is, however, con- 
venient to introduce the new quantities 


v=¥ “i 7 
v'=Ws—vs J 


YW’ and w’’ are the streamfunctions for the thermal flow 
in the lavers between levels 1 and 3 and levels 3 and 5, 
respectively, and are of course measures of the mean 
temperatures in these lavers 

Applving finite differences and using the boundary 
conditions for w, we arrive after certain simple lanipula- 


tions at the following set of equations for w: and ay: 


, 2f,2 {2 / eas ) 
oV*w, WoT p? 4 P V? V “Vy ) 


Pp? 
. V Ve’ V’-Vn 


Ww 


Vive 


} 20 





Va — w . w V2(V...-Vy’’ 
}/? /? P ¥ 
V ‘we V’''’.Vn EW’. VE" 
P is a constant pressure interval equal to 3314 eb Equa- 
tions (2.6) is the svstem which has to be solved for @, : d 


», given the three streamfunctions Ys, ¥’, and y’’ 
The three streamfunections will be prescribed by the 


CXPPessions: 


y Uy + Ay sin ka 
y’ ('y A’ Sih hea + @’ ee 
y’’ U’’y+A” sin (kr+a’’ 


U;, U’, and U” 


thermal winds in the lavers above and below this level 


are the zonal wind at level 3 and the 


All three will be assumed to be constant. Ay, A’, and 21” 
are the amplitudes of the streamfunctions. It is seen 
that Ay, kA’ and kA” are the maximum meridional wind 
components in the three fields. a’ and a’’ measure thie 
phase lag of the two thermal fields relative to thie 
streamfunction at level 3. Positive values of the ) 
mean that the thermal field is lagging behind the stre: 
field at level 3. 

It is pertinent to mention that the values of a’ and 
determine the mean slope of the trough and ridge lines ni 
the lavers between the levels 1 and 3 and the levels 3 and 


5, respectively. a’ and a’’ determine together a measure 


of the change in slope with height of the pressure systen 5 
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The flow pattern at level 1 is obtained by adding the 
expressions for y; and y’. 
formulas we may write the streamfunction at level 1 in 
the form: 


Vv; (03+ U7 )y+ A, sin (kr +7;) (2.8) 
where 
A, (A,7+ A”’?+2A,A’ cos a’)! (2.9) 
and 
tan y= Ama. (2.10) 
Az T A’ cos a’ 


It is thus seen that y, will be positive if a’ >0 provided 


ew is not so large that A,+ A’ cos a’<—0. In a similar 


wav the flow pattern at level 5 can be written: 


y { U’’)y+ A; sin (kr—y; (2.11) 
where 
| A;* . 2.4,A”’ cos a’ 2.12 
id 
A’’ sin a’ 
tan y ;* 2.13) 


A;—A’’ cos a’ 


(gain we find that Y and a” in general have the same 


sign for atmospheric flow patterns. We have thus seen 


that positive values of a’ and a’’ mean that the flow at 
level 1 is lagging behind the flow pattern at level 3, while 
the flow pattern at level 5 precedes that at level 3. Posi- 


tive values of a’ and a’’ indicate therefore the usual west- 
ward tilt of the pressure systems, but the slope of a system 
is not necessarily linear in pressure. 

{ 


vel the following set of 


When the expressions 2.7) are substituted in the rig! 


hand sides of equations 2.6) we 


equatlons 


~ 2] - / oA key ‘) 
Ww p? Wo-T p2* P és h 4 COS Ad 
hA’(8—L’k*) eos (kaa 
> (2.14 
c yA oP : / 1713-4. eos ka 
pe pe pP = adie - 
kA’’(B+U'’k?) cos (ka+a’’ 





ll 


The solutions to the system (2.14) will be of the form: 


w RB COS her B, COS hea t a’) B COs (hea a’’ 
715 
(', cos kr+-C, cos (kr+a’)+C, cos (kr+ al’ J 
he amplitudes B and (5, can be determined by 
substitution of (2.15) into (2.14) and equating the coeffi- 


When 


find the following 


eMlents of cos kr. COs (h 4 a’), and COS (ker t a’’). 
th procedure is carried through we 


expressions for the two vertical velocities: 
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Applying simple trigonometric 
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The quantity, A, appearing in (2.16) and »17) isa 
notation for the following expression 
», oh Te if i.\" : 
a—(o+ 28) ose 2B)—(fY. cats 


The meridional components 7,=h.ly cos ka, 0! = kA’ 


COs 


hey to’), and v’’ kA’’ cos hea a’’) have also been miro 
duced in the solutions (2.16) and (2.17) 
From the knowledge OL ws and Ws obtained from (2.16 


and (2.17) we ean evaluate the divergence aut the levels l. 
3, and 5, if we at the same time use the boundary condi- 
tions for w; 1.e., w=0 for p=0 and p=p 

from the 


We vel 


The divergences are evaluated continuity 


equation in finite difference form. 
v-V) iF 


AND SOME 


DISCUSSION OF THE SOLUTION 
EXAMPLES 
the 


vertical motion zonal cross-sections are appropriate. To 


To best illustrate distribution of divergence and 
produce such cross-sections it 1s hecessary lo obtain values 
for the parameters, which are parts of the coefficients 

In the first example we have chosen a wavelength of 
1000 km. The stability 
determined in such a way that they correspond to the 


It has been 


parameters, o, and o,, were 
usual increase of this parameter with height. 


(Wiin-Nielsen [6]) that the 


with 


with 
the 


shown variation of o 


pressure is described cood accuracy using 


expression 








TABLE 1 Values of Ll’ and Ll” 
Dat GMI n i 1 
Fet S, 1050, 1200 22. | 14. 
Fet 10, 1959, 1200 21.4 12.6 
Fet 12. 1959, 1200 21.2 1S 
Feb. 16, 1959, 1200 17.2 08 
Feb, 20, 1959, 1200 xs 13.2 
Me Is 12 
50 \? 
» 
o\p) a; ( ) (3.1) 
p 


where o, 1s the value of the stability at 50cb. The present 
operational model applies a value of ¢; equal to approxi- 
mately 3 MTS-units. Using (3.1), we get 

a9 MTS-units, oyx2.25 MTS-units. (3.2) 


The the thermal winds, 0” and 7’, were 


determined from actual wind data for a number of indi- 


values of 
vidual days. The mean zonal winds were available at the 


levels 85, 50, and 30 eb. ©’ and 1’ were then obtained by 


linear extrapolation and interpolation. The values are 
given in table 1. 
The mean values, (’=18 msec. U"’=12 m.sec.~ 


were used in the examples. The striking feature for the 


five days in February 1959 is that the average shear is 
50) percent larger in the upper laver than in the lower 
laver It is interesting to see whether this result holds 
in general for the middle latitude in winter. From wind 


stutisties for the period February 10, 1959 to April 10, 


1959 it was found that the mean value of 0” was 
16 m.see.~' and of U7’ 9.5 m.see. Data presented by 
Petterssen ({4], p. 98) for Larkhill covering a 2-vear 
period gnve the values [ 14 m.sec.~!, U’’=7 msec. 
It may therefore be safely concluded that the vertical 
shear of the horizontal wind is somewhat larger above 
than below 50 cb. in middle latitudes in winter. The 
maximum values of the meridional wind components were 
chosen to be: 7’=12 m.see.~', v’’=8 m.sec and pr 
16 msec 


The last for which numerical values 


have to be chosen are a’ 


two parameters 
In the first computa- 
were selected. These 


14° as computed 


and a’’. 
tion the values a’=10° and a’’=20 
and Y 


5360 


values correspond to y,=4 
(2.10) (2.13). As 
km. we find that the trough line at level 1 is only about 


from and corresponds to 4000 
14 kin. behind the trough line at level 3, while the trough 


line at the lower level (level 5) precedes the trough at 


level 8 by about 156 km. This first case has therefore 
a very small slope but a somewhat larger slope in the 
lower layer than in the upper. This difference in slope 
is Characteristic for many atmospheric systems. 

Figure 2 shows the distribution of the divergence in a 
zonal cross-section through a half wavelength. The heavy 


dashed line in the central part of the figure is the position 


The heavy solid line is the isoline for 


of the trough line 
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1000! Lo 
cross-section 
1000 
km. and the stability parameters correspond to the usual increase 
of stability with height. West is to the left, east to the right. The 
heavy dashed line is the position of the trough line. 


FiguRE 2 Distribution of divergence in a zonal 


through a half wavelength. In this case the wavelength is 


this 


We that the 


divergent level outside the trough region is situated a 


zero-divergence. find in case non- 


little above the 500-mb. surface. However, as we ap- 
proach the trough line from the west the line for zero- 
divergence goes down to lower levels and intersects level 
5 to the west of the trough position. Approaching the 
trough line from the east the zero-divergence line goes up 
and intersects with level 1 to the east of the trough at 
this level. We find therefore convergence at all levels in 
the trough. At the lower level we find convergence to 
the east of the trough and divergence to the west, while the 
opposite is the case at the higher level. 


The distribution of divergence shown in figure 2 


agrees 
to a very large extent with the distributions computed 
by Charney [1], Fleagle [2], and Hinkelmann [3]. It 
seems therefore that the major features of the distribution 
of divergence in a continuous model can be reproduced 
by a three-parameter model. 

It was found of interest to investigate the distribution 
of divergence in the same zonal cross-section, if the vertical 
variation of the stability parameter was disregarded. — In 
this computation it was assumed that a where 


The resulting dis- 


04 0 
og; Was set equal to a standard value. 
tribution of divergence is reproduced in figure 3. Com- 
paring figure 2 with figure 3 one sees that the distributions 
are similar in a qualitative sense. The divergence pattern 
at level 3 has, however, an amplitude about 6 times 
larger in figure 3 than in figure 2, while the divergence 
is increased by a factor of 2 at level 1. It seems therefore 
important at least to incorporate an increase of the sta- 
bilitv parameter o with height in an operational model 
because the divergence and the vertical velocity are quite 
sensitive to this variation. 

This fact is also illustrated in figure 4 which shows the 
distribution of the vertical velocity through a half wave- 
length at level 2. The solid line corresponds to the case 
where the vertical variation of stability is incorporate l, 
while the dashed line is the distribution of the vertical 
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FIGURE 3. 
section as in figure 2. In this case the stability parameter is set 


crTross- 


zonal 


Distribution of divergence in the same 


to a standard value and does not vary in the vertical. 


velocity in the case where o,.=o4=03. The damping in- 
fluence of the increase of the stability parameter, ¢, is 


clearly seen in the figure. A similar figure (not repro- 


duced) at the lower level (level 4) shows a very small 
difference between the corresponding curves. This is 


explained by the very small vertical variation of o in the 


lowest laver of the troposphere. 


4. ON SOME FACTORS INFLUENCING THE 
MID-TROPOSPHERIC DIVERGENCE 


The divergence at level 3 estimated from the values of 
w, and w,, which in turn are found by solving the two 
coupled w-equations, will be used in a numerical prediction 
in the divergence term in the vorticity equation. It is 
therefore of importance to investigate the different factors 
which influence the distribution and magnitude of this 
divergence. In the present model with three information 
levels we have at least a first approximation to the curva- 
We have 


further found in the preceding section that the vertical 


ture of the vertical profile of the zonal wind. 


variation of the stability parameter is of importance for 
the magnitude of the mid-tropospheric divergence, but 
not so much for the distribution. 

In order to simplify the first part of the discussion we 
shall for a moment neglect the vertical variation of the 
stability o. With this simplification we obtain from (2.16) 


and (2.17) 


f 
Ue ) : ‘ | A bhaa. 
¥ of? +-3f 7/k? a(l 


(U’—p/k*)r’ —(U + B/k*)v"" (4.1) 


The first term in the bracket of (4.1) represents the in- 
fluence on the divergence patterns at level 3 caused by the 
crude measure of the curvature of the vertical profile of 
It is seen that (07 l 


the horizontal wind. ’’) is propor- 
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FIGURE 4 Distribution of vertical velocity through a half wave- 


level 2 (3345 eb.) 


case with vertical variation of stability, 


The solid line corresponds to the 
the dashed line to the 


length at 


case of no vertical variation of stability “T” designates position 


of the trough line 


tional to the second derivative of Ll’ with respect to pres- 


sure: 
i Oe | ' 
~— (U’—U"’). (4.2) 
dp? p? 
If re Sia aus in the example represented in figures 2 
and 3 the first term will result in convergence where 


r, >0, 1.e. between the trough and the next ridge down- 
stream, and divergence between the ridge and the following 
trough. habe 
The influence of the first term is therefore to increase the 

pe 


’’ but the term has no 


The opposite distribution will result if 0” 


instantaneous speed of propagation of the wave if 0” 
and to decrease the speed if (’< 1 
influence on the instantaneous deepening or filling 

In barotropic forecasts for 500 mb. it is often found 
that a trough is forecast to move too slowly in the southern 
portion, essentially in the region of the subtropical jet 
According to synoptic investigations of the ver- 


find U’>U"” 
It is therefore possible that the contribution from the first 


stream. 
tical wind profile in this wind system we 


term in (4.1) will help to remove this error 
Let us next turn our attention to the last two terms in 
(4.1). Bik?) 1s 


positive we find that both terms produce convergence 


If the wavelength is so short that (1 


are positive and divergence where they 
are the 
Y’ and y’’, are lagging behind the flow 
terms contribute to produce convergence in the trough 


where p’ and p’’ 


negative. If therefore temperature patterns, 


¥;, both of these 
and thus give a tendency for deepening. However, for 
B/k*) 


second term changes sign and the two last terms counteract 


waves which are so long that (0 is negative, the 


each other. In a qualitative wa) this argument shows 
that the shorter waves develop faster than the longer 
Equation (4.1) is somewhat simplified because we have 
ussumed that O20, If this assumption is removed we 
get a more complicated expression for the mid-tropospheric 
2.16) and (2.17) we obtain 


divergence. Subtracting 
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fk? ; - — _ 
(V°V) i. U’ (ak? + f,2/P?) —U "" (ook? + fo?/P?) } 1 
(o,k?+ f,7/P?) (UL —B/k* Dv’ 
(a.k?4- f,2/P2)(U'’ + B/k?)v’"|. (4.8) 


The more general expression (4.3) shows that it is not 
only the difference in the vertical shear which counts in 
the first The thermal wind €’” in 
the upper laver is now multiplied by a factor depending 
is multiplied 
by a factor containing the stability in the upper layer. 


term in the bracket 
on the stability in the lower laver, while U7’’ 
Using stundard values of 0, and o, we Can find the ratio 


| ted ie? 
first term in the bracket of (4.3). 


which would give no contribution from the 
The ratio (U’/U’’). is 
given as a function of wavelength in table 2 computed 


from the formula: 


(7? ok? f2/P? 
( z) : ry (4.4) 
[ o,k*-4 Jo P 

the table it that the 
of stability has the greatest influence on the divergence 
With no variation of the stability 


From is seen vertical variation 


for the shorter waves. 
we get convergence ahead of the trough if only (U’’—U'’’) 


that 07’ 


have to be at least two times larger than U’’’ to get the 


is positive. The values in table 2 show would 
Saine sign of the divergence for a wave with a wavelength 
of 4000 km. 


modification is 


The table shows also that only a minor 
the 
derived 


introduced for long waves as 
the 


For a given value of (U’/L’ 


very 


compared to results from the simplified 


formula (4.1). ’) character- 
istic of atmospheric conditions, say 1.5, the first term in 
(4.3) would result in divergence ahead of the trough in 
short waves, but convergence in the same region for long 
Wives. 

With respect to the last two terms in (4.3) we find some 
modification of the magnitude, but no change of the sign. 
The contributions from the terms are, however, large 
enough to give convergence to the east of the trough 
line in figure 2 as well as in figure 3. 
the 
zonal 


Figure 5 shows distribution of convergence and 


divergence in a cross-section In a more extreme 


case. The distribution shown in figure 5 was computed 
using the same parameters as in figure 2 except that the 
phase differences between the temperature fields and the 
stream function at 500 mb. were larger in the case il- 
lustrated by figure 5. In the construction of figure 5 it 
that a’ $5 and a’’ 9) 
that the temperature field in the upper laver lags 's of a 


was assumed which means 


TABLE 2. Values of (U' U''). as a function of wavelength. 
L&<1W0-* km l 2 $ i i i) 7 sS te 12 i4 2s 
u'/( $f 0 | 2.4 | 2.0 iva l 1.4 1.3/1.2 Lae) oe 1.0 
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Figure 5 Distribution of divergence in the same zonal cross- 


section as in figure 2. In this case the phase differences between 


the temperature fields and streamfunction at 500 mb. were larger 


than in case shown in figure 2. 


wavelength behind the streamfunction at 500 mb., while 
The 
result of the assumption is a larger slope of the trough 
From the formulas (2.10) and (2.13) it follows that 


the lag in the lower layer is 4 of a wavelength. 


line. 
the trough at level 5 precedes the 500-mb. trough by 
about 300 km., while the trough at level 1 lags about 
210 km. behind the 500-mb. trough. 

The main effect of the greater vertical tilt of the system 
is, as seen from figure 5, to produce a larger slope of the 
The 


of the trough line has now a maximum in the midtrop- 


convergence pattern. convergence in and ahead 


osphere amounting to about 2 107° sec.~' with divergence 
above and below. 


5. REMARKS ON THE DIVERGENCE IN LONG 


WAVES 


VERY 


Due to the special difficulties encountered in fore- 
casting the very long waves in the atmosphere it is of 
interest to investigate the distribution of divergence in 
such waves in the present model. It is apparent from 
(4.3) that 


become of minor importance when the wave number is 


the terms containing the vertical stability 
small. In the same case it is seen that the two 8-terms 
become dominating in the last two terms of the bracket in 
(4.3). 


sien and therefore tend to compensate if 7’ and v’’ are of 


These two 8-terms appear, however, with opposite 
the same order of magnitude and in phase. The impor- 
term for the determination of the di- 
vergence distribution in_ this the difference, 
(C’—U"’). As we find that (0”’— WU’) is positive in the 


jet stream regions where the waves have the greatest 


tance of the first 
model is 


amplitude, the contribution is to produce convergence to 
the east of the trough and divergence to the west in the 
very long waves and thus obtain a decrease of the retro- 
gression. Forecast experiments with the model will show 


whether the magnitude of the divergence is large enough 


to control the very long waves. 
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We may estimate the divergence from (4.3). Assuming 
the same values of the stability as in figure 3 (o,=9 and 
o,=2.25 MTS-units) and further (’=18 m.see.-', U’’ 
12 m.sec.~'! we find in middle latitudes (f)=107* see.~') 


lor ke? 0.2 > 
the hemisphere: 


10-"m.~*, corresponding to two waves around 


v-V); 10.29 v3—3.78 v’+6.29 v’’|K107*> see.-! (5.1) 


Assuming further that the waves 73, 7’, and v”’ are in 


phase and adopting the values r3,max= 16 m.see.~', 76 ..=12 
sec.~', and vf4,=S8 m.sec.~! we find that 
(V°V).~—107' cos kr (5.2) 
if the streamfuncetion at level 3 is 
¥3;= —U; y+ As sin kv. (5.3) 


We find therefore in this case a distribution of divergence 
which would counteract the retrogression of the very long 


waves. 
6. SUMMARY AND CONCLUSIONS 


The distribution of vertical velocity and divergence in 
a three-parameter model of the atmosphere has been 
investigated. Section 2 contains the formulation of the 
problem and the formal solution of the equation for the 
vertical velocity in the model for simple sinusoidal waves. 
The solution obtained in section 2 is discussed in the next 
section, where the importance of the vertical variation of 
static stability and the vertical profile of the horizontal 
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wind is pointed out. Examples showing the distribution 
of divergence in zonal cross-sections are shown. It is 
found that the distribution of divergence in the three- 
level model to a large extent resembles distributions by 
others with more complicated models. 

The different factors influencing the mid-tropospheric 
It is found that 
the vertical profile of the wind is most important for the 


divergence are investigated in section 4. 


speed of propagation of the waves, while the vertical 
slope of the waves contributes to the development 

The last section contains a discussion of the distribution 
of divergence in very long waves and an example is 
computed. 
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ABSTRACT 


An equation for vertical velocity 


velocity and divergence in special cases characterized by pronounced failures of barotropic forecasts. 


is applied to data from 850, 500, and 200 mb. for calculations of vertical 


The results 


of the calculations show that a non-development situation was adequately described by the equivalent-barotropic 


picture of a single quasi-horizontal surface of non-divergence. 


The onset of mid-tropospheric development, as shown 


by the appearance of large errors of the barotropic forecast, was characterized by the appearance of a double sur- 


face of non-divergence, with a deep mid-tropospheric convergence layer in the vicinity of the trough line 


This 


picture of development is confirmed by a second case study 


The appearance of the mid-tropospheric convergence layer is related to the low-level cold push, the tilt of the 


flow patterns with height, and the high-level jet stream, confirming synoptic studies by J. J 


others 


1. INTRODUCTION 


In considering the problem of baroclinic development, 
usually associated with cyclogenesis, one is often led to 
consider the various theoretical studies of baroclinic in- 
Such studies for prediction 


stability. two-parameter 


models vield quantitative instability criteria. However, 
when applied to atmospheric data, two-parameter predic - 
tion models produce forecast failures strikingly similar to 
those of barotropic forecasts, especially in cases of pro- 
nounced baroclinic development. 

In considering the capabilities of a two-parameter 
model which contains representations of vertical shear 
and vertical velocity, from which vertical momentum ad- 
vection can be computed, it is difficult to imagine one 
which does not contain, by unplication at least, a quasi- 
horizontal surface of non-divergence. It therefore seems 
likely that a feature of the atmospheric circulations not 
found in two-parameter models, geostrophic, balanced, or 
primitive equation models, night be of considerable im- 
portance in the development process. One such feature 
which should be considered is the mid-tropospheric di- 
vergence, which can be computed from an atmospheric 
model having three or more parameters in the vertical. 

This study will attempt to throw some light on the 
development) process by means of computing vertical 
veloc ities and Varlous derived fields by means of a three- 
level model of the atmosphere. For the purpose of this 
study, develoy ment can be considered as the processes 
which lead to error in the barotropic forecast (over areas 
of good data We shall also exclude special mountain 
and frietion effects from this study. In order to bring 


out more clearly the development processes, the actual 


events will be deseribed in terms of the errors of the baro- 





George, H. Riehl, and 


tropic forecasts. Vertical motion, divergence, and ver- 
tical momentum advection can then be discussed in terms 


of these errors. 
2. THE VERTICAL VELOCITY EQUATION 


The vertical velocity equation is Obtained in the usual 
method from the vorticity equation and the adiabatic 
equation written in the following form: 


oc Ow 
——t- Ye ° 0 (1) 
ae V-Tn—n op 
0 (Od Oo 
+ . WwW 0). (2 
ai ap TVV ap to ) 


In (1) and (2) V is the horizontal wind; ¢ is the relative 


vorticity ; n is the absolute vorticity; nm is the absolute 
VOrtcity at 500 mb. which can also be used to represent 
a Vertical average of n; w dp dt, the vertical velocity ; 
is the geopotential; and o is a measure of the static sta- 
bility given by o adlné/op, where @ is specific volume 
ind @ is potential temperature. 

We now assume that for the purposes of calculating 
horizontal advection the horizontal wind is non-divergent, 
and can be streamfunction, 1. @., 


V=kxVy. 


strophic approximation necessary to obtain an equation 


represented by a 


In place of the usual formulation of the geo- 
for w, we shall make an approximation to the temperature 


field that 0¢/ 0p=/foy/Oop. 


mations into (1) and (2) leads to the w-equation: 


(fo ; ov] 
) [= (V-Vn) V2 (V-0 - )] 


This involves the assumption that horizemtal variations 


Introduction of these approxi- 


iny Ow 


V-w+ rs 
o Op” 
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of static stability can be neglected, although ¢ is permitted 
to vary in the vertical. The appearance of fm as a 
coclficient of the second term in equation (3) was required 
for special reasons not relevant to this study. However, 
different versions of equation (3) were used in calculations 
in which (a) static stability varied freely, involving extra 
terms in o, and (b) f; (at 45° N.) replaced fm. The results 
which applied to this study were not significantly affected 
by these changes. 

The finite difference system in the vertical consisted of 
using data from 850, 500, and 200 mb. An interpolation 
was made enabling the use of 800, 500, and 200 mb. in the 
calculations. In computing finite differences at the low 
levels, the variable height of the ground was taken into 
The upper and lower boundary conditions 
V,-Vp, at the lower 
boundary p,, the standard atmosphere pressure at the 


consideration. 
consisted of w=0 at p=0 and w 
variable height of the ground. 

For finite differences in the horizontal the Joint Numeri- 
cal Weather Prediction Unit’s (JNWP) grid on a polar 
stereographic map having a mesh length of 381 km. at 
60° N. In the computation of Laplacians 
and Jacobians, a consistent system was used in which the 


was used. 
components wu and ¢ of the wind were taken over a double 
mesh length. The calculations were made on the IBM 
704 computer. 

Before proceeding to the case studies, let us consider 
certain implications of the vertical velocity equation. 
For this purpose we can simplify equation (3) to the form: 


O*w 


ap ea (4) 


where A The 
right-hand side of the equation is represented as the 
forcing function, F’, determined by the vertical variation 
of vorticity advection and the Laplacian of the tempera- 


k?o/ fn, k being a horizontal wave number. 


ture advection. In a quasi-geostrophic model, n-1 values 
of F can be obtained from n parameters in the verticai. 
The usual two-parameter prediction model therefore con- 
tains sufficient information for only one value of F in the 
If we then take F as independent of pressure, 
0 at 1000 mb. and 
at O mb., the w-profile is symmetric about 500 mb., where 
This is represented 


vertical. 
and apply as boundary conditions @ 


a surface of non-divergence is found. 
by figure la. Of course it would be possible, in the 
absence of other than climatological information, to make 
an assertion regarding the variation of F with pressure in 
such a model. The only effect of this would be to change 
the pressure of the non-divergent surface. However, a 
continuous non-divergent surface corresponding with an 
isobaric surface would still exist. 

Suppose instead that we consider a model with sufficient 
parameters in the vertical (minimum of three) to provide 
information on the slope of the vertical F-profile. Figures 
lb and le represent the shapes of the w-profile associated 
will two simple distributions of F. It is clear that if the 
slope of the F-profile changes sign from one place to 
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Ficure 1.—-Three vertical profiles, (a), (b), and (c), of forcing 
function, F, and vertical velocity, w, as a function of pressure p 
PND refers to the pressure at which no horizontal divergence 


exists. 


another significant areas of mid-tropospheric convergence 
or divergence may occur. An atmospheric model with 
three or more levels is thus suitable for calculation of mid- 
tropospheric divergence. Although detailed 
calculations will be necessary before the question of the 


numerical 


number of levels required for adequate vertical resolution 
can be settled, it is clear that three levels represents the 
minimum number required if mid-tropospheric divergence 
is an important baroclinic phenomenon. 


3. CASE STUDIES 


(I). JANUARY 21-22, 1959 


On January 22, 1959, a notable cyclogenesis was 
observed over North America. 
the Great Lakes and travelled 
deepening rapidly. Of special interest was the fact that 
large errors rather suddenly appeared in even short-range 


A sea level cyclone passed 


across into Labrador, 


barotropic forecasts. 

The development was preceded by the motion of a 
500-mb. trough through the central United States. Figure 
2a shows that 24 hours in advance of the development the 
atmosphere was behaving barotropically, as shown by the 
small errors of the barotropic forecasts. The cross-section 
of horizontal divergence for the same time (fig. 2b) shows 
the familiar equivalent barotropic pattern of a quasi- 
horizontal surface of non-divergence near 500 mb. with 
compensating convergence and divergence patterns above 
and below the surface. Some distortion of the pattern is 
found over the mountains as a consequence of the ascent 
of the air to the west and the descent to the east of the 
highest. terrain. 

Figure 3 shows a more complicated situation by 1200 
GMT of January 21. This could be regarded as a period of 
transition from an equivalent barotropic situation to a 
baroclinic situation. The errors of the barotropic forecast 
for the subsequent 12 hours increased to amounts which 
were not negligible. The cross-section (fig. 3b) shows a 
split of the surface of non-divergence in the vicinity of the 
trough line. The signs of the convergence and divergence 
patterns in mid-troposphere are consistent with the errors 
observed in the barot ropic forecast. 
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Fiaure 2.—(a) Error in deeafeet (dashed lines) of the 12-hr. barotropic forecast from 0000 Gut, January 21, 1959, superimposed on the 
500-mb. contours (solid lines) for that time. (b) Cross-section of horizontal divergence in units of 10~-® see.-'. The numbers in 


circles are w in units of 10-5 mb. see. Heavy dashed line is the trough line. 























Figure 3.—Same as figure 2 except for 1200 Gut, January 21, 1959. 
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Figure 4.—Same as figure 2 except for 0000 Gut, January 22, 1959. 


In figure 4 the development is shown proceeding in 
The error of the barotropic forecast over the 
eastern United 
fig. 4b) shows a most significant change from conditions 
In the vicinity of the trough line there 


strength. 
States was extreme. The cross-section 
24 hours earlier. 
was a double surface of non-divergence, one at high levels 


and one at low levels, with a deep layer of convergence in 


mid-troposphere. This cross-section of a developing 
trough resembles similar cross-sections obtained — by 


Charney [3] and by Hinkelmann [6] from idealized data, 
as well as those computed from real data by Bundgaard 
[2] and by Fleagle [4]. 

Having the vertical velocities, one can compute the 
effects at 500 mb. of the vertical advection of momentum. 
The tendency of the 500-mb. height arising from the ver- 
tical momentum advection is shown in figure 5. It can 
be seen that in this situation this effect is generally of a 
sign. to reduce 500-mb. error, but can account at the most 
for about 20 percent of the error. 

In considering that vertical velocities of 10 
or more over large areas are relatively common in winter, 
that vertical displacements of 100 mb./day 
Should result. If the vertical velocity is downward in a 
trough line, this could result in a strengthening of the flow 


* mb./see, 


one sees 








Contribution of the vertical advection of momentum 
units of ft. per 12 


Figure 5. 
to the 500-mb. height tendency (dashed lines; 
hr.) for 0000 amr January 22, 1959, superimposed on the 500-mb 


contours (solid lines) for that time 
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through a deep layer by 20 percent in 24 hours, an amount 
which is by no means negligible. 

In figure 6, the intersections of the isentropic surfaces 
with the cross-section are shown. These show the pro- 
nounced cold dome, and cold push, to be a feature of the 
lower two-thirds of the atmosphere. Above this we find 


a Warm push, associated with the low warm stratosphere. 
(II). DECEMBER 6, 1959 


On December 6, 1959, a cyclogenesis occurred along the 
eastern United States coast. Figures 7 and 8 show the 
extreme error of the barotropic forecast over the south- 
eastern States as the low-level cyclogenesis took place. 
The computed vertical velocities at 650 aad 350 mb., 
together with the computed 500-mb. divergence, are shown 
in figures 9, 10, and 11. Examining the areas of the 
southeastern States, we can see that a large area of 500-mb. 
convergence Was associated with two separate features. 
Ahead of the trough the ascending motion Was stronger at 
high levels than at low levels, and in the trough line the 
A 500-mb. 
divergence area is found just behind the trough line, where 
The 500-mb. 
divergence pattern of figure 10 can be compared with that 


descending motion was stronger at low levels. 
the descent was strongest at high levels. 


of figure 6b in the paper by Brown and Neilon [1], who 
obtained their map by assuming that the entire error of the 
barotropic forecast was a result of the 500-mb. divergence 
pattern. Although there are considerable differences in 
details, there is a good resemblance of the major features. 

Another view of this can be obtained from cross- 
sections running just north of the main centers of activity. 


The orientation of the cross-section (fig. 12) was chosen 
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Cross-section for OOOO Gm, January 22, 1959, showing potential temperature in the same vertical plane as in figure 4b 


to coincide with that presented in the following paper 
by Brown and Neilon. This gives the reader an oppor- 
tunity to check two completely independent determina- 
tions of vertical velocity and divergence. The cross- 
section of divergence (fig. 12b) shows a pattern in the 
developing trough similar to that of figure 4. The vertical 
velocity cross-section (fig. 12a) shows that in this case the 
convergence associated with the trough was associated 
mainly with a strong ascending motion at high levels. 
This relatively strong high-level ascent, as well as the 
strong descent behind the trough, was a consequence of 
a strong 200-mb. jet stream in strongly curved flow 
through the trough. If there had been no phase shifts in 
the vertical, this would merely have resulted in a faster 
vastward displacement of the 500-mb. trough than given 
by the barotropic forecast, since at 500 mb. there would 
have been convergence ahead of and divergence behind 
the trough. However, the phase shift in the vertical 
resulted in an extension of the 500-mb. convergence area 
into the trough line as well as an enlargement of the 
convergence area. This conclusion is in agreement with 
the results of Wiin-Nielsen [9] who showed that the lagging 
of the temperature field behind the flow leads to a mid- 
tropospheric tendency for deepening in the trough line 


4. GENERAL REMARKS 


A feature of both cases was the tendency for the sinking 
of cold air to be at a maximum in the lower levels. A 
well-known characteristic of strong cold-air outbreaks is 
the tendency for the maximum cold-air advection to 
occur in the lowest half of the atmosphere, with a low, 
prominent in the higher levels. 


warm stratosphere 
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Figure 7.—Error in deeafeet of 12-hr. barotropic forecast from 
0000 Gmr December 6, 1959. 
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Ficure 8.—Error in decafeet of 12-hr. barotropic forecast from 
1200 amr December 6, 1959. 


According to the analysis of the w-equation presented 
earlier, such cold-air outbreaks therefore ought to be 
characterized regularly by a very strong low-level di- 
vergence, a relatively low surface of non-divergence, and a 
laver of mid- and upper-tropospheric convergence. In 
such cases, failures of barotropic forecasts are often 
observed. Cyclogenesis is also frequently seen under such 
conditions. In support of this diagnosis we can cite the 
evclogenesis study by George and collaborators [5], who 
found the strong low-level cold-air push to be a reliable 
precursor of cyclogenesis. They found the 850-mb. chart 
to be most useful in recognition of the low-level cold 
push. J. Austin* has pointed also to the occurrence of 
the abnormally low warm stratosphere during the initial 
stuges of cyclogenesis. 

* In lectures at two American Meteorological Society meetings at Washington, D.C. in 


195} and 1955, 











Ficure 9.—850-mb. contours (solid lines) and 650-mb. # (dashed 


lines) in units of 10 


> mb. see. 


! for 1200 cmt, December 6, 1959 
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Fiaure 10.—500-mb. contours (solid lines) and divergence (dashed 


lines) in units of 10°° see 


! for 1200 Gut, December 6, 1959 
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Figure 11.—200-mb. contours (solid lines) and 350-mb. w (dashed 


lines) in units of 10 


3 mb. see. 


' for 1200 cmt, December 6, 1959 
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Figure 12.—Cross-section (a) vertical velocity (10-% mb. see.~!) 


Wiin-Nielsen [9| showed in his 3-level analysis of a model 
disturbance that an upward increase in vertical wind 
shear leads to a pattern of 500-mb. divergence which 
tends to speed up systems (convergence ahead of a trough, 
divergence behind), as compared with barotropic motion. 
Taba’s |S] typical vertical profiles of the subtropical jet 
stream show the normal existence of such conditions. 
It is a common experience at the JNWP Unit that the 
barotropic forecasts move features too slowly in the 
areas occupied by a subtropical jet stream. If a sub- 
tropical jet stream is found over the forward side of a 
trough, as in the schematic representation of figure 13, 
the resulting mid-tropospheric convergence accentuates 
the development. Here we should point out that Riehl 


(7! observed this association in 1947. His paper contains 





and (b) divergence (10~° see 


T 
25°N 59°W 


!), for 1200 amr, December 6, 1959 


special mention of the upward increase of shear under 
the high-level jet stream, in association with the occurrence 


of cyclogenesis. 
5. CONCLUSIONS 


The equivalent barotropic representation of the atmos- 
phere, with a quasi-horizontal surface of non-divergence 
near 500 mb. separating divergence patterns of opposite 
sign above and below, was a realistic representation of 
actual conditions during a period preceding cyclogenesis. 
However, during two different cyclogenetic situations, 
the equivalent barotropic idea became invalid. During 
these development periods a double non-divergent surface 
in the of the developing troug!. 


appeared vicinity 
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Fiagure 13. 


Between upper and lower divergence regions a deep mid- 
tropospheric convergence layer appeared. The invalida- 
tion of the equivalent barotropic representation was 
further demonstrated by the appearance of large errors 
in the barotropic forecasts. 

The particular atmospheric features responsible for 
this mid-tropospheric convergence area near the develop- 
ing trough line were the following: 

(1) The low-atmospheric strong cold-air advection in 
and behind the low-level trough surmounted by high-level 
warm air advection. This led to a vertically asymmetric 
w-profile, with low-level sinking and mid-tropospheric 
convergence. 

2) The normal phase shift with height of the system, 
in which the sinking cold air is brought in at lower levels 
under the upper cyclonic system. This accentuates the 
effects mentioned under (1) above. 

3) The participation of a strong high-tropospheric jet 
Stream in the strongly curved upper flow. This accentu- 
ales the pre-trough convergence in mid-troposphere. 


Schematic representation of factors associated with mid-tropospheric development 


These features are represented in figure 13, in which the 
low-level cold push, the high-level jet stream, and the 
phase shift with height are indicated. The inset w-pro- 
files represent conditions at the locations A and B. One 
cannot expect that all of these contributing factors will 
be present in each case of mid-tropospheric development, 
but the errors of the barotropic forecasts will reach the 
largest values when all are present in strength. 

The most serious limitations on the accuracy of the 
calculations presented here are the use of the geostrophic 
approximation, and the lack of a better vertical resolution 
of the atmosphere. It is just in development situations 
where the vertical velocities and ageostrophic motions 
reach their greatest magnitudes. At the same time the 
smaller-scale horizontal motions become of greater im- 
portance to the developmental process. It is therefore, 
a very large step from diagnosis to prognosis. However, 
in deriving an w-equation of the general type used in this 
study it is necessary to state some kind of wind law. In 
view of the well-known sensitivity of numerical predictions 


82 


to the wind law used, it seems that a generalization of the 
wind law to be used in the w-equation to that of the balance 
equation might be of sufficient interest to warrant recod- 
ing the machine calculations. 

The problem of vertical resolution of the calculations 
is aggravated by the presence of a tropopause of variable 
height. If the vertical wind shears above and below the 
tropopause are each relatively invariant (but of opposite 
sign) serious errors can be introduced into the estimates 
from streamfunctions at 


wind derivatives 


widely separated levels on opposite sides of the tropo- 


of vertical 


pause. Further experiments to study this question seem 


to be in order. 
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ABSTRACT 


With the aid of an electronic computer, case studies of wind analyses at the 850-mb., 700-mb., 500-mb., 400-mb., 


300-mb., and 200-mb. pressure levels have been made. 


from the « and v wind-component analyses are investigated. 


The divergent and non-divergent wind components resulting 


For the cases considered, the streamfunction fields 


are slightly superior to the Joint Numerical Weather Prediction operational fields, obtained initially through use of 


the “‘balance equation.”’ 


from the winds by previous investigators employing hand-analysis techniques. 


The magnitude of the horizontal wind divergence values are comparable to those obtained 


However, the divergence patterns 


are not sufficiently accurate for the strict requirements necessary for numerical weather forecasting 


1. INTRODUCTION 


Analyses of several meteorological parameters have been 
made during the past several vears at the Joint Numerical 
Weather Prediction Unit (JNWP) through use of elec- 
tronic computer techniques. At the present time analyses 
of the heights of the 850-mb. and 500-mb. pressure sur- 
faces are made twice daily on an operational, routine 
basis.’ This analysis scheme was described recently by 
Cressman [1]. 

Many attempts have been made in the past to obtain 
accurate representations of the flow at constant pressure 
surfaces, but in most cases these attempts have been made 
Results of 
studies made in recent vears such as those by Landers [4], 
Murakami [5], Rex [7], and Taba [11] point to the pos- 
sibility of obtaining direct, accurate wind analyses over 
regions of dense data coverage. These authors point to 
considerable skill in obtaining the horizontal wind diver- 
gence at several levels in the atmosphere for their indi- 
vidual case studies. It is a well-known fact that 
with sufficient data coverage the usefulness of divergence 


through use of hand-analysis techniques. 


even 


computed from instantaneous wind observations is de- 
pendent upon a high degree of observing accuracy. 
Although one can hope that observational errors are 
random and thereby can be reduced considerably by 
modern analysis procedures, it is questionable whether 
the reduction in error will be enough to produce sufficiently 
accurate divergence patterns to be useful for future fore- 
Casting purposes. 

The purpose of this study was to determine whether the 
wind fields at constant pressure surfaces in the troposphere 
and lower stratosphere could be accurately analyzed over 
the INWP octagonal grid in an objective sense through 
use of a sufficiently large electronic computer. 

\n attempt was also made to determine to what extent, 
il any, the geostrophic assumption, made in the JNWP 


——————_——__ 


December routine basis the 


temperature 


5, 1960, JINWP began analyzing on an operationa 


ind wind fields for 850, 700, 500, and 300 mb 
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height analysis program, is damaging to the wind forecasts 
from the operational divergent one-parameter forecast 
model. 


2. ANALYSIS PROCEDURE 


All computations were performed by the IBM 704 
electronic computer on all or part of the regularly ar- 
ranged JNWP octagonal grid which covers most of the 
Northern Hemisphere north of 10° latitude. The grid- 
point interval was 381 km. at 60° latitude on a polar 
stereographic projection. The analysis scheme which was 
employed was a modification of the JNWP operational 
height-analysis routine which is described in detail in [1]. 
Only the major features of this program will be mentioned 
here. 

This analysis technique utilizes an initial estimate of 
the pattern of the field being analyzed and adjusts this 
guess in a prescribed manner to fit the data. Therefore, 
for meteorological data which are sparse in certain large 
regions of the grid, it is desirable to use a first guess which 
closely approximates the final analysis over the entire 
region. 

Attempts were made in the early stages of this study 
to obtain first estimates of the wind fields at all levels 
under consideration, excluding 500 mb., through use of 
a linear regression equation employing the 500-mb. stream- 
function wind, obtained from the JNWP operational “ bal- 
(9],. and the 500-850-mb. geostrophic 
The resulting extrapolated wind fields 


ance equation” 
thermal wind. 
contained large influences of the 500-mb. cvclostrophic 
winds which were improperly positioned. This was pur- 
ticularly noticeable below the 500-mb. level where the 
tilt with height of the pressure patterns was more pro- 
nounced. 

The non-divergent part of the geostrophic wind (for 
further details of this see [8]) of the operational height 
analyses was tested for its usefulness as a first guess 
The resulting analyses were not sufficiently removed fron 
because of the 


the geostrophic influences. Therefore, 


sensitivity of the final non-divergent and, particularly, 
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the divergent wind components to an accurate first wind- 
field guess, it was decided to use the balance-equation 
winds for the first estimates of the final analyses at all 
levels. 

All wind data were checked in the JNWP operational 
automatic data Four 
passes are made through the fields, and during each of 
these scans the following correction ( to the wind com- 


processing system.’ successive 


ponent wis made at each grid point: 
C=—WE. 


Here, / is the error of the linearly interpolated u-wind, 
obtained from the first guess (or from the values com- 
puted in the previous scan), at the observation location, 


and 
N2—@? 


Wate 


Here, d is the distance from the observation to the grid 
point being modified, and N is the radius of a circle cen- 
tered at the grid point. N varies from 4.75 grid intervals 
during the first sean to 1.0 grid intervals during the fourth 
sean. W is set equal to unity during the final sean. If 
more than one observation falls within the circle pre- 
scribed by NV, a correction C is computed for each obser- 
vation and the average correction is applied at the grid 
point under consideration. 

In addition a five-point (center and four immediate 
surrounding grid points) smoothing operator is applied 


after scans 2 and 3: 


= a 
U= Up V7u- (1) 
T; 


Here, u denotes the smoothed value, uw is the value at 
the central point, and y? is the single grid-increment 
finite difference Laplacian operator. 

After the last sean is completed a weak 9-point smoother, 
described by Shuman [10], is applied. 

After the uw wind component was analyzed by the above 
method, the procedure was repeated for the ¢ wind com- 
ponent. 

Because of the large barotropic 500-mb. height errors 
and the existence of dense upper-air data coverage Over 
the region of interest, the cases for 0000 Gur, November 
16, 1959 and 1200 amr, December 6, 
Winds were analyzed over North America at the following 
pressure levels: 400, 300 and 200 mb. In addition, the 
850-, 700-, and 500-mb. levels were analyzed over the 
entire JNWP octagonal grid. The scale factor of the map 
was used throughout this study. 
to the earth rather than to the map were obtained and 
The effects of 
considered in this study. 


3. SYNOPTIC SITUATION 


Since the results of the tests of both cases studied were 


1959 were studied. 


Thus velocities relative 


examined. truncation errors were not 


A vertical consistency check is made on all winds greater than 15 kt. at or near the 
mandatory reporting levels from 1000 to 200 mb. A detailed description of this test has 


been given by Dent [3] 
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similar, it was decided to direct most of the discussion (o 
the 1200 emt, December 6, 1959 case. 

At this time the synoptic situation was as follows: a 
surface Low existed over Lake Huron with a cold front 
extending southward through western Florida and then 
southwestward through the Gulf of Mexico (fig. 1). A 
warm front extended eastward from this Low and then 
northeastward along the St. Lawrence Seaway. A weaker 
cyclonic circulation was centered about 450 miles east of 
Delaware. High pressure existed over the western United 
States with a ridge extending to a second high cell centered 
over southern Texas. A strong cyclonic circulation which 
had moved eastward through western Canada was cen- 
tered near Ft. Nelson, British Columbia. 

By 0000 Gut, December 7, 1959, the Low over Lake 
Huron had weakened into a trough which extended to a 
new cyclonic circulation centered over southeastern Vir- 
ginia. Twelve hours later this storm was located in 
northeastern Pennsylvania and had intensified consider- 
ably—central pressure was slightly less than 980 mb. 
By this time the storm in western Canada had moved into 
the western edge of Hudson Bay and increased slightly 
in circulation. Aloft a strong ridge existed over the 
western United States with a sharp, intense trough over 
the Mississippi River valley at 1200 Gur, December 6, 
1959 (fig. 2). This trough moved rapidly eastward during 
the following 24 hours. 


4. EVALUATION OF RESULTS 


According to Helmholtz’s theorem the horizontal wind 
vector (WV) may be separated into a rotational non- 
divergent part and an irrotational divergent part; L.e., 


V=kxVy+ Vx. 


Here, k is the unit vector directed upward, V is the hori- 
zontal gradient operator, and yY and x are the stream- 
function and velocity potential, respectively. It was 
decided to test each of these wind components separately. 


A. THE NON-DIVERGENT WIND FIELDS 


From the analyzed uw and ¢ wind components it was 
possible to obtain the relative vorticity (¢) from the finite 
difference form of the following equation : 
Goo, Ww OU 
fo Or Oy 


7 


In equation (2) g is the gravitational acceleration, fo 
the Coriolis parameter at 45° latitude, and 


Thus, ¥ has dimension length. The winds were differ 
enced over a double grid increment. 


From equation (2) it was possible to obtain the field of 


the streamfunction, given the boundary values. The-e 
values were obtained by assigning the streamfunction of 


an arbitrary boundary point the value of the height of t! « 
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Fiagure 1.—850-mb. divergence (solid lines) and surface isobars 
dashed lines) 1200 Gort, 1959. D 
and C refer to divergence and convergence centers, respectively. 

units of 10-7 see.-! indicates 


and fronts for December 6, 


Values are in Cross-hatching 


precipitation areas. 


pressure surface obtained from the JNWP height analysis. 
From the wind components normal to the boundary it was 
possible to obtain a first estimate of the streamfunction at 
the next boundary point. Thus estimates were obtained 
foreach boundary point. Since the integrated divergence 
over the region was not necessarily zero, a large dis- 
crepancy between the y-gradient between the first and 
last point and the analyzed winds at those points could 
Therefore these boundary p-values were adjusted, 
This 


result. 
and the integrated divergence was also obtained. 
was used later in obtaining the velocity potential. 

The sealed streamfunction field was then obtained by 
the Liebmann relaxation technique using a double grid- 
increment finite difference Laplacian.’ The 500-mb. 
result for 1200 cur December 6, 1959 is shown in figure 2. 
The streamfunction isolines appear to fit the wind data 
quite nicely. The isotachs were computed from the 
analyzed winds which contained both the divergent and 
non-divergent components, rather than the streamfune- 
tion winds, to illustrate the accuracy of the analysis and 
at the same time to illustrate the effects of the analysis 
For the 0000 cmt, November 16, 
1959 ease, the smoothing routines after scans 2 and 3 of 
Therefore, this 
produced analyses which were in closer agreement with 
the data. 

The 500-mb. non-divergent wind fields obtained from 


smoothing routines. 


the analysis procedure were omitted. 


this procedure were then verified over the United States 


laxation using a single grid-increment finite difference Laplacian produced 


imilar rel 
i streamfunction values which differed from the double grid-increment ones by 


s 50 feet at 500 mb. in the vicinity of deep Lows 
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500-mb. streamfunction lines and observed winds for 
1200 aut, December 6, 1959 


FIGURE 2. 


rhis 


non-divergent wind field was obtained from the stream- 


and southern Canada using the reported winds. 


function field through use of centered finite differences 
over a double grid-increment. The winds at the grid 
points were then linearly interpolated to the observation 
locations and compared to the observed data. The pro- 
cedure was also applied to the JNWP operational stream- 
function, from the The 


root-mean-square (RMS) and average vector errors are 


obtained balance equation. 


noted in line 2 of table 1. The same verification pro- 
cedure was applied to the total analyzed wind (line 1, 
table 1). 
smoother used in the December case. 

It is not surprising that the RMS vector errors of the 


Here note the effects of the 2d and 3d sean 


non-divergent wind are less than the operational errors. 
This is due to the fact that the operational streamfunction 
winds are obtained from a height analysis which utilizes 


TABLE 1.—Rool-mean-square wind and height verification  sta- 
fistics 
O000 owt. Nov. 16, 1959 mw omeTr De 6, 1050 
500-mb. fields 
INWP Wind INWP Wind 
1. Analyzed winds (Knots 5.3 4.5 09 7.8 


Initial non divergent 
winds (knots 15.3 °(12.2 11.4 (9.3 13.8 (11.3 12.0 1s 


to 


3. 12-hr. barotropic forecast 

winds (Knots 15.9 (18.3 14.8 11.8 2.3 In. 0 1 f 
4. 24-hr. barotropic forecast 

winds (Knots 2.2 (6.9 19.3 16.2 $0.3 (25.3 27 22 
5. 36-hr. barotropic forecast 

winds (knots 22.9 (19.7 22.2 (18.8 $2. 27.8 20.8 (2 
6. 36-hr. y-forecast vs. y 

verifying winds (Knots 31.7 20.3 $1.2 28.8 
7. 36-hr. g-wind difference 

knots 13.0 5.1 
&. Initial heights (feet 44 SI 71 Tt 
). 36-hr. barotropic fore cast 

heights (feet $45 $05 123 43 

*Values in parentheses are the aver tor er 








approximation to fit the grid-point 


In addition, these winds 


the 
heights to the observed winds. 


geostrophic 


were used for the first guess in obtaining the wind analyses 
made in this study. 

The 500-mb. streamfunction field the 
analyzed winds was further tested by making a barotropic 
The JINWP operational 


The resulting 


obtained from 


forecast from this initial field 
forecast program was used for this purpose. 
forecast streamfunction winds were compared with those 
obtained from the routine barotropic forecasts issued by 
INWP by the same verification procedure described above. 
(See lines 3, 4, and 5 of table 1.) Figures in line 6 are the 
36-hr. RMS vector errors computed from the 36-hr. fore- 
cast streams and the INWP verifving streams. <A 36-hr. 
computed the JINWP 
Notice that in all forecasts the 


persistence measurement from 
streams is given in line 7. 
barotropic forecast winds made from the wind analyses 
are a slight improvement over the JNWP operational 
results. 

The height fields were obtained from the streamfune- 
tion fields in the usual manner from the balance equation. 
After linear interpolation of these values to the observa- 
tions, the RMS height errors, computed over the same 
regions as were the winds, were obtained and are presented 
in lines Sand 9 of table 1. The results of the larger initial 
height errors from the wind analy ses were surprising since 
it was expected that the geostropic assumption made in 
the INWP height analyses was damaging. If the balance 
equation and its method of solution were a complete repre- 
sentation of the atmosphere, a streamfunction which 
exactly represented the observed winds should produce a 
field which exactly the 

Although the balance equation produces a fric- 
tionless, the streamfunction 


values at the boundary points are set equal to the height 


height represented observed 
heights. 
non-divergent wind and 
values, it is difficult to explain this diserepaney quantita- 
tively through these approximations inherent in the bal- 
ance equation and its necessary boundary restrictions. 
Phillips [6], while studying the Appalachian storm = of 
November 1950, noted a similar discrepancy at 400 mb. 
The 36-hr. forecast-minus-verifving height pattern from 
the INWP forecasts for the December case is presented 
in figure 3. A similar pattern was obtained from the 


wind-analvsis forecasts. 


B xENT WIND FIELDS 


THE DIVEK 

The fields of horizontal wind divergence were obtained 
directly fromthe wand ¢wind analyses at all of the levels un- 
derconsideration through use of the centered finite difference 
approximation. The result obtained at 850 mb. for the 
December case is presented in figure 1. The region of 
convergence Newfoundland 
ward over the St. Lawrence Valley is reasonable since con- 


found over und southwest- 


siderable amounts of precipitation were observed in this 
region. Convergence was also found over the southeast- 
ern United States where evelonic development occurred 
this level over 


in the succeeding hours. The results at 


the western portion of the continent are unreliable because 
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Figure 3.-- 500-mb. JNWP operational 36-hr. forecast-minus-veri- 
fving height pattern valid at 0000 Gut, December 8, 1959. Values 
are in units of 10° ft 
: ; . s ; 

of the presence of the mountains. The area of con- 


vergence near the southern tip of Texas is questionable. 

The velocity potential field was obtained from the 
divergence field in the same manner that the streamfune- 
tion field was obtained from the relative vorticity field; 


that is, by inverting 


itiect— 
hy Or OV 
where 
. Jo 
x x 
q 


is the sealed velocity potential. A Liebmann relaxation 
was performed on the 500-mb. data again using a double 
Rather than 


assign specific boundary values for X it was decided to use 


grid-increment finite difference Laplacian. 


the average boundary normal gradient of x which was 
aVailable from the procedure described in A above. 

The result of this computation for the December case 
The corresponding 
patterns of The 
divergent wind is directed perpendicular to the isolines 


for 500 mb. is presented in figure 4a. 
divergence are shown in figure 4b. 
from centers of divergence (positive values) to centers 
of convergence (negative values). Over a large portion 
of the central United States this wind is less than 
Winds between 5 kt. and 10 kt. are found in the western 
and eastern United States and in southern parts of the 
These latter features may be due, in part, to the 
Noti 


the extremely large divergent wind speeds which resulte:| 


5 kt 
grid. 
effects of the first guess used in the wind analysis. 


over the south-central region of Canada. These resu 
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from the large-amplitude, short wavelength, divergence 
paitern obtained from the wind analysis. 


C. STREAMFUNCTION TENDENCY 
AND VERTICAL VELOCITY COMPUTATIO 


The vorticity equation can be written in the following 
consistent, simplified form (see Wiin-Nielsen {12}): 


oy . fi 
' —J (pn) —— VV. (3) 
v ot M n q ) 
Here J is the Jacobian and 7 is absolute vorticity. The 
500-mb. sealed 12-hr. streamfunction tendency _ field 


resulting from the individual terms in the right side of 
the above equation and their sum are presented in figure 
the The latter 
fig. 5d) is an approximation obtained from the balance 
12 hr. 
baro- 


5 together with “observed” tendeney. 
equation streamfunction fields 12 hr. before and 
1200 «mr, 1959. The 12-hr. 


forecast made from the wind analysis study 


after December 6, 


tropic and 
the operational streamfunction > showed positive height 
errors over United States. 
from figure 5a, the contribution of the first term of equa- 


the eastern As can be seen 
tion (3) produced negative stream tendencies over this 


region which were too small in magnitude. The magni- 
tude and pattern of the streamfunction tendency resulting 
from the divergence term (fig. 5b) demonstrates that if 
these divergence fields were to be included in a numerical 
fields directly to 
the 


initial tendencies would undoubtedly contain large errors 


forecasting model which used these 


btain initial streamfunection tendencies, resulting 


over regions Where equation (3) was valid. 
At this point it 
modified 


that another 


was 


should be mentioned 
the 

streamfunction 
the relative 
over a large region is conserved in the following equation: 


form of equation used to 
estimate the 


Wiin-Nielsen has also shown that 


vorticity 


local tendency patterns. 


vort icity 


a 
,O fy fon 
gM boy, ry VY, (4) 
Or g q 

V, and 7 refer to the total wind and absolute vorticity, 
The tendeney field 
puted from this equation was very similar to that found 


respectively. streamfunction com- 


using equation (3). However, the contributions of the 


lividual terms could not be evaluated by the relaxation 


procedure because of large ‘“‘pillows” of opposite sign 
which resulted from each of the forcing functions. 
This result may be more clearly understood by con- 


sidering the equation 


: a) he 
2 di — Ve ALA, 
ae 1 73," V, 


JA . 


Where Ais the Theoretically for a sufficiently 
larve area the above integral approaches zero provided the 
absolute vorticity and the diverge nee are uncorrelated. The 
field the 
divergence term of equation (4) contained large positive 


500 mb. 


area 


streamfunetion tendency which resulted from 


errors at This was evident from the discrepancy 
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500-mb. velocity 


wind speed 


1200 December 6, 1959. (a 


(solid lines 


Fieure 4 GMT, 


potential and corresponding divergent 


isotachs (dashed lines). Wind speeds are in knots. Plus and 
minus refer to centers of divergence and convergence, respective ly 
b) Field of horizontal wind divergence from which (a) was 
obtained Units are 10°-7 see Dashed line indicates line along 


which vertical cross section of figure 7 was taken 


found in the 


which was streamfunction tendency fields 

resulting from the divergence terms of equations (3) and 

(4). Since the errors of a were positive over the central 
ot a 

regions of the grid, YW were negative. 


the errors in V° 
or 


Thus, according to the above equation, the vorticity and 








A / 


7. 
it Jiv.n b)— v-V: and (¢ J(¥,n)- v-V: (d) The streamfunetion tendenev obtained from operational streams 12 hr. - 
/ q 


( 


fore and 12 br. after 1200 Gut, December 6, 1959 Units are 10° ft./12 hr 


divergence obtained from the wind analyses were posi- 
tively correlated In other words, at 500 mb. in regions 
ol evelonic VOrUcity horizontal divergence was prevalent, 
and the reverse was generally true for regions of anti- 
evelonic vorticity 

At the 850-mb. level the individual contributions of the 


terms of equation (4 produced tendency pillows of the 
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igure & The contribution to the 500-mb. 12-hr. sealed streamfunction tendencies for 1200 cmv, December 6, 1959, fr 








opposite sign to those found at 500 mb. Thus at these 
levels the vorticity and divergence were negatively cor 
related. Similar results were found at 700 mb. 

These general results are in agreement with our preset 
knowledge of the atmosphere. Assuming that the ab- 
solute vorticity patterns are similar throughout the lower 
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half of the troposphere, the results imply that the level o! 
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PaBLE 2.— Mean absolute values of divergence (10-7 sec. —') 


Pressure level Nov. 16, 1959 Dec. 6, 1959 


(mb 0000 GMT 1200 GMT 
S50 36 43 
700 35 36 
500 38 17 
400 46 64 
300 is 63 
200) 43 59 


least divergence over the octagonal grid exists somewhere 
between the 700- and 500-mb. pressure levels. This is 
further substantiated by the averaged divergence values 
computed over the United States at each of the pressure 
levels. (See table 2.) 

A further test was applied to the divergence patterns 
which resulted from the wind analyses at 500 mb. If the 
long-wave stabilizing term and the mountain term, which 
exist in the JNWP operational barotropic forecast model, 
are neglected, the barotropic forecast equation reduces to 


A 


2 


Dt —J (We,ne). (5) 


The subscript B refers to the scaled barotropic stream- 
function. Subtracting equation (5) from (3) produces the 


following approximation: 


A 
Here y, is the scaled barotropic streamfunction forecast 
error (forecast minus observed) and Af is the time incre- 
ment. The average wind divergence field obtained from 
the 24-hr. JNWP barotropic forecast from 0000 Gar, 
December 6, 1959 is presented in figure 6. This diver- 
gence pattern approximates the 24-hr. average pattern 
necessary to correct the operational JNWP barotropic 
forecast. The general features of this pattern may be 
compared with those of figure 4b, obtained directly from 
the wind analyses. 

The approximated vertical velocity patterns were ob- 
tained from the divergence patterns (smoothed) of the 
analyzed winds through the use of the continuity equation: 

Ow 
Op 
Here p refers to pressure and w=dp/dt. 
w is negligible at 1000 mb. and that the divergences vary 
linearly between the levels at which they were measured, 


—VeV. 


Assuming that 


it is possible to compute w at all of the data levels. A 
cross-section of these patterns for the December case is 
presented in figure 7a. Figure 7b is the corresponding 
cross-section of the divergence fields. These cross-sections 
were made along a line extending from Norfolk to Los 
Angeles (dashed line in fig. 4b). These 


sunilar to those of the classical picture; i.e., convergence 


patterns are 
at ‘ower levels and divergence at the upper levels appear 
in regions ahead of the trough and the reverse in regions 
to the rear of the trough. The vertical velocity errors 
introduced by the neglect of mountain effects are un- 
Cross-sections of vertical velocity 


do ibtedly significant. 
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L : 
FiGcuRE 6.—500-mhb. estimate of the horizontal wind divergence for 
1200 cmt, December 6, 1959, obtained from the error field of the 
INWP 24-hr. operational streamfunction forecast from 0000 G7, 
December 6, 1959. 10-6 see.~! 





Units are 


and divergence fields from Cressman’s [2] three-parameter 
quasi-geostrophic model for this time and over the same 
region may be used for comparison. (See fig. 12 of Cress- 
man’s paper.) The large discrepancy found in the vertical 
velocities near the Continental Divide is undoubtedly due 
to the lower boundary restriction used in the present 
studv. However, the large-scale features of both sets of 
cross-sections are in agreement in the main. The magni- 
tudes of the divergence values are generally larger in the 
wind analysis study. 


5. SUMMARY AND CONCLUSIONS 


Although the w and ¢ wind fields were analyzed inde- 
pendently, the wind fields produced by the analysis scheme 
“fitted” the wind observations more accurately than did 
the original wind fields that were used as first estimates 
Since it was desirable to have a 

the final analysis, the JNWP 


operational streamfunction, obtained from the balance 


of the final analysis. 
good first estimate of 


equation, was used for this purpose. 
The the 
analyses also fitted the data more accurately than did 


streamfunction winds obtained from wind 


the JNWP streamfunction winds, and in both cases the 
former produced slightly improved barotropic wind fore- 
casts. However, the initial height field produced from 
the wind-adjusted streamfunction field, which in turn 
was obtained through use of the balance equation program, 
was inferior to the operational JNWP height analysis. 
The employment of the newer streamfunction in pref- 
the JNWP operational one is not 


erence to present 
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economically feasible with the present JINWP numerical 
that a 


the same quality as obtained 


forecasting models. Preliminary tests indicate 


streamfunetion of about 
in this study may be produced from the balance equation 
streamfunction by a more economical procedure. 

Briefly, this entails a slight modification of the present 
machine analysis program to enable it to produce stream- 
function analyses, using wind data only, through use 
of the following equations (rather than the geostrophic 


assumption) : 


oy oy 
u — ’ r . 
OV Ou 
This improved streamfunetion and the one found in the 
present study Tian also prove to be more desirable, parti- 
cularly in the upper troposphere and lower stratosphere, 


since the ellipticity criterion (¢, 7-5? where Cg IS the 


g 


geostrophic relative vorticity), which must be satisfied in 
solving the balance equation, That be violated in regions 
of dense data coverage 

The average magnitudes of the horizontal divergences 
resulting from the « and ¢ wind analyses are comparable 
to those found by Landers and others using hand-analvsis 
techniques. The magnitudes obtained by direct measure- 
ments appear to be larger than those obtained by most 
implicit calculations. As a first approximation, if we 
assume that the 500-mb. barotropic forecast errors are due 
largely to the horizontal wind divergence contribution, it 
would appear that the amplitudes of the divergence fields 
obtained in this studs are of the correct order of magni- 
tude. Although difficult to evaluate, the divergence pat- 


terns over flat) terrain uppear to be more accurately 
represented at the lower levels of the troposphere where 
observational density and accuracy are greater. 
However, the divergence patterns obtained in this study 
could not be utilized by a numerical forecasting model 
Which ineluded the evelogenetic mechanism and required 


accurate initial horizontal wind divergence measurements. 


Cross-section obtained from the observed winds of 1200 Gur, December 6, 1959, 
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A NOTE ON CLEAR AIR TURBULENCE DURING APRIL AND MAY 1960 


FRANK J. SMIGIELSKI 


t Meteorological ¢ 


[Manuscript received July ¢ 


1. INTRODUCTION 


During the months of April and May 1960 a number 
of teletypewriter transmissions of pilot reports of clear 
air turbulence in the United States were gathered for a 
The study was undertaken 
to see if any clues to forecasting clear air turbulence in 
the levels 500 and 300 mb. 
These levels were chosen because of the larger number 
of flights that are made at levels between 18,000 and 
The 


approach of this study was to outline areas where clear 


study of this phenomenon. 


between could be found. 


30,000 ft. both by jet and turbo-prop aircraft. 


air turbulence was encountered by pilots and to examine 
these locations in relation to upper-air flow patterns, 
horizontal and vertical wind shears, and 500-mb. height 


changes. 


2. CLEAR AIR TURBULENCE 
IN RELATION TO UPPER-AIR FEATURES 


Clear air turbulence was most often reported when the 
500-mb. well 
In this type of circula- 
tion, the turbulence was quite often reported to the east 


flow showed a closed Low or a defined 


trough as its primary circulation. 


and southeast of the trough line or low center and im- 
mediately to the west of the area of maximum winds and 
Figure | 
shows a typical 500-mb. pattern when clear air turbulence 
Wis reported, In 


Was associated with strong horizontal shear. 


clear air turbulence was 
17,000 21,000 ft.; all 
reports indicated moderate to severe turbulence. Strong 


this 


cause 


reported at levels between and 
horizontal and vertical shear was present, as is shown in 
figure 2, a cross section from Peoria, Il. to New York 
City based upon upper-air observations 2 or 3 hours 
before the time of the turbulence reports. 

Figure 3 shows a case of closed circulation at 300 mb. 
that was associated with reports of clear air turbulence. 
Turbulence was reported at 29,000 ft. in the vicinity 
of Morgantown, W. Va. about 2 hours before the time 
of the upper air observations. Earlier in the day, severe 
clear air turbulence was reported at 31,000 ft. between 


Cleveland and Chicago. Figure 4 shows another type 


of upper-air pattern that gave rise to clear air turbulence. 

Balzer and Harrison [1] noted a lack of cases of clear 
al turbulence over the Rockies and northern Plains. <A 
possible explanation of this could be the mean positions 


Mice, U.S 


Weather Bureau, Wa 


hington, D. 


1960: revised Dec. 21, 1960] 


of troughs and ridges aloft over the United States: a 
mean ridge lies over the Rockies, and a mean trough 
farther east over the Central States [2]. This 
give a greater amount of clear air turbulence over the 


would 


area from the eastern Plains eastward to the east coast. 
Along the Pacifie coast, it is not uncommon in certain 
portions of the year to have a cut-off Low over the south- 
western United States, and along the California coastline, 
with a trough northward along the coast. In these cases, 


a greater amount of clear air turbulence cases would 


occur in the western United States. Clem [3] showed a 
clearcut maximum level of clear air turbulence over the 
southwestern and northwestern United States between 
27,000 and 40,000 ft. in winter, and this would fit the 
pressure patterns for clear air turbulence. 

In the study of horizontal wind shear in relation to 
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Figure 1.—500-mb. chart, 1200 Gut, April 5, 1960. Solid lines are 


contours: dashed lines are isotherms; solid circle indicates report 


of clear air turbulence. 
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Figure 2.—Vertical section from Peoria, Ill. to New York City, 
1200 Gur, April 5, 1960. Solid lines are isotachs; dashed lines are 
isentropes; blocks indicate reports of turbulence; circle indicates 


report of smooth air. 


clear air turbulence, it would seem that the direction 
of travel of the aircraft would be a factor in the effect 
of the turbulence on the aircraft. In those reports of 
turbulence for which the direction of aircraft travel could 
be determined it was found that the direction was most 
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Figure 3.—300-mb. chart, 0000 amr, May 9, 1960. Solid lines are 
contours; solid circle indicates report of clear air turbulence. 
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Figure 4.—500-mb. chart, 1200 aur, April 8, 1960. Solid lines are 
contours; dashed lines are isotherms: solid circle indicates report 
of turbulence. 
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TEMPERATURE LAPSE RATE, °C PER 5000 FT 
Figure 5.—Clear air turbulence reports during April and May 1960, 


plotted as a function of vertical wind shear and temperature lapse 


rate, 


often from the Jow pressure to the high pressure side of 
and most often was nearly per- 
An added effect in cases 


the upper-level flow 
pendicular to the streamlines. 
where a well-developed trough existed, especially at 500 
mb., was upslope motion, as suggested in figures 1 and 
+ by the isotherms. 

Arakawa [4] derived an equation that vields a critical 
value of horizontal shear of 15 knots per degree of latitude. 
Any value above this limit would vield heavy turbulence. 
Harrison and staff [5] give a critical shear value of 50 
knots per 150 nautical miles. It is thought, however, 
that beeause clear air turbulence is generally reported in a 
narrower horizontal band than 150 miles, the practical 
calculation of shear from the wind difference over a finite 
normal distance should perhaps use an interval of about 50 
nautical miles. This would yield more representative 
values of the shear necessary to produce clear air turbu- 
It is unfortunate that the wide spacing of upper-air 
reporting stations makes use of this smaller interval 


lence, 


extremely difficult. 

In an attempt to correlate vertical wind shear and 
vertical temperature gradient over a 5,000-ft. interval in 
the turbulence cases, a great variation in the paired values 
was found. Figure 5 shows pilot reports of clear air 
turbulence plotted with ordinate of vertical wind shear 
over a 5,000-ft. interval regardless of sign and abscissa of 
change in temperature in degrees Ce!sius over the same 


vertical interval. The cases studied make up only a small 
suniple, but they do indicate that clear air turbulence can 
exist with a vertical wind shear of less than 10 knots per 
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5,000 ft., and with temperature lapse rates of 5° C. per 
5,000 ft. or greater. The graph shows that a vertical wind 
shear of 15 knots per 5,000 ft. is enough to produce clear 
air turbulence if the lapse rate is more than 5° C. per 
5,000 ft. The strong vertical wind shear of 55 knots per 
5,000 ft. that is plotted on the graph (fig. 5) is for the same 
date and time as figures 1 and 2. 

Clear air turbulence was reported in areas of 24-hour 
height falls at 500 mb. 
little over 0.5, but suggests a relation of kinetie height 


The correlation here was only a 
changes to clear air turbulence. Taking the changes over 
a 24-hour period tends to mask such a relation; it is 
believed that if changes over a smaller interval of time, say 
6 hours, could be computed, a better relation between clear 
air turbulence and height changes at 500 mb. might be 
found. 


3. SUMMARY 


Clear air turbulence was most often reported at levels 
above 15,000 ft. whenever the 500-mb. chart showed a 
closed Low or a well defined trough. The turbulence was 
then found in the southwesterly flow ahead of the trough 
line, on the low pressure side of the maximum wind field. 

The reports received indicated the existence of clear air 
turbulence in the vicinity of the jet core, on the poleward 
or low pressure side of the jet; this is in close agreement 
with the findings of Bannon [6] and Jones [7] in investiga- 
But clear 
air turbulence is not a phenomenon of the jet stream alone, 


tions of clear air turbulence in the British Isles. 


for it was reported as far as 20,000 ft. below the level of 
maximum winds. 

There was evidence that a strong horizontal shear across 
a trough line adds to the intensity of clear air turbulence. 
Some evidence pointed to the existence of clear air turbu- 
lence in the right quadrant of negative height falls at 500 
mb., but the correlation was not very strong. 
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ABSTRACT 


Clear air turbulence (C-A-T) occurrences over the United States reported by pilots in April 1960 are analyzed 
with reference to weather conditions, wind speed, horizontal and vertical wind shear, thermal stability, and jet stream 
location and curvature. The day-by-day plot of the occurrences shows that the peak of the C-A-T activity was 
associated with sharp and well developed troughs which extended far to the south. In the analysis of the data, it is 
found useful to separate the cases east and west of 103° W. longitude. The cases west of 103° W. indicate greater 
importance of higher stability and less importance of large horizontal wind shear. The cases east of 103° W., away 
from the influence of the western mountains, indicate greater importance of decreased stability when combined with 
either large horizontal or vertical wind shear. Well developed jet streams were apparent with the occurrences in 
both areas. However, the western cases show a preponderance of cyclonic curvature while the eastern cases show 


nearly equal division between cyclonic and antieyelonie curvature. 


1. INTRODUCTION TABLE 1.—Clear air turbulence occurrences over the United States for 


April 1960. 


The purpose of this study was to make a preliminary 
Light 195 


analysis of the Clear Air Turbulence (C-A-T) occurrences Sendarane 3 
for April 1960 with reference to the meteorological con- Talpabus tikeisaine "6 


ditions in order to determine the problems involved in 


obtaining all the necessary data, in analyzing the data, 
and in interpreting the results. ‘This study provides some In this analysis, the parameters were evaluated as near 
rroundwork and guidelines for a broader investigation of — as possible at the actual location and level of the reported 
('-A-T in addition to furnishing some tentative conclusions | turbulence even though greater values of these parameters 
regarding the meteorological conditions associated with might be found at some nearby location or level. The 
('-A-T actual evaluations were based on the upper-level charts 

Pilot reports for this period collected from the Wash- as analyzed at the National Weather Analysis Center 
ington National Airport, United Airlines, and American (NAWAC) and on all the appropriate rawinsonde, 
Airlines, were seanned for clear air turbulence (C-A-T) radiosonde, and pilot balloon reports. In most cases, 
occurrences. Table 1 is a tabulation of all cases of | this analysis called for considerable interpolation between 
observation sites and times. As far as possible, a spacing 
of 100 miles centered at the site of actual occurrence was 
used in the evaluation of the horizontal wind shear and 
a 2000-ft. layer centered at the level of the actual 
2. METHOD OF ANALYSIS turbulence was used for the evaluation of the vertical 
wind shear and the stability. 

Many of the reports were not specific as to whether the 


turbulence over 10,000 feet, except those which definitely 
were described as in or near thunderstorms, squall lines, 
or with icing in clouds. 


Since light turbulence is not a major concern to aircraft 
operations, only the moderate and severe turbulence , ; on 
' . _ : turbulence occurred in clear air or in clouds. The 
occurrences were used in the analysis. The following , ' 
ear : , presence of thunderstorms, squall lines, widespread 
parameters were examined in connection with each ° , 
showers and cloud cover as determined by surface maps 

occurrence ¢ ? > al 
——— , and other weather data reduced the number of C-A-1T 
a. Surface weather—thunderstorms, squall lines, ; 4 
cases to approximately 170. 
fronts, showers, clouds. ; 


b. W ind Speed. . 3. SYNOPTIC PATTERNS 

c. Horizontal wind shear (in knots per 100 miles). 

d. Vertical wind shear (in knots per 1000 feet). The locations of the C-A-T occurrences were plotted 

e. Thermal stability (potential temperature lapse on both the maximum wind and 300-mb. charts for thie 
rate, ° C. per 1000 feet). chart time nearest the time of the turbulence. Forty 


f. Location, intensity, and curvature of the jet percent of the C-A-T occurrences were between 0600 and 
stream. 1800 GMT and 60 percent between 1800 and 0600 Gm. 
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The larger number of occurrences in the latter period is 
likely due to a greater number of flights during these 
afternoon and evening hours. 

figure 1 shows the day-by-day occurrences of C-A-T 
during this period. The four major peaks in the C-A-T 
activity corresponded to the presence of large-amplitude 
major troughs aloft in the C-A-T areas during these 
periods. This is in agreement with the findings of 
Harrison [1]. 

This is illustrated in the sequence of 300-mb. charts as 
analyzed by NAWAC for the April 12-14 period (fig. 2). 
As the trough deepened and moved eastward, the C-A-T 
occurrences moved with the trough until it began to 
dissipate. Similar sequences were apparent during the 
other peak periods but are not shown in this paper. 


4. LOCAL PARAMETERS 


The first examination showed a wide scatter in the 
values of the various local meteorological parameters in 
connection with the observed C-A-T occurrences. It was 
not until the occurrences over the coastal ranges and the 
Rocky Mountains were separated from the occurrences to 
the east of the Rockies that some patterns and conclusions 
could be drawn. 

Table 2 shows the range in the values of the wind speed, 
horizontal and vertical wind shear, stability, and the 
presence and curvature of the jet stream associated with 
the C-A-T occurrences. A line at approximately 103° W. 
was used to separate the cases over the western mountains 
and those to the east of the mountains. 

The scatter in the values of the wind speed and the 
vertical wind shear is about equally great in the eastern 
and the western occurrences. Over 70 percent of the 
occurrences had wind speeds of 50 knots or more and about 
45 percent of the occurrences were associated with vertical 
wind shear of 4 knots per 1000 feet or more. <A great 


ranite 2.— Number of cases of occurrence of values of meteorological 


param ters associated with the C-A-T occurrences. 


Wind speed (knots Horizontal wind shear (knots per 


100 miles) 
O49 Mw-74 75-09 100 o-4 15-24 25-34 > 35 


East of 103° W IS 28 19 s Is 24 20 12 
West of 103° W 27 33 17 i 36 34 15 5 


Vertical wind shear (knots per |) Lapse rate (potential temper 


1000 feet) ature) (°C. per 1000 feet) 
0 1.1 2 10 7.0 0 a 2 40 6 
1.0 2.4 3.9 6.9 1.0 2.4 3.4 9 
East of 103° W 13 6 13 13 20 26 12 Is 1! 3 
West of 103° W 21 10 9 21 15 20 11 15 Is 12 


Jet stream present Jet stream curvature 
Cold Low 
present 
Yes No Cyclonic Anti 


eyclonic 


K f 103° W 63 10 34 29 3 
W of 103° W 76 14 a) 20 6 
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Day-by-day occurrences of clear air turbulence 


(C-A-T) during April 1960 


FiGure | 


preponderance of the occurrences were associated with a 
pronounced jet stream. However, the western cases 
show a preponderance of cyclonic curvature while the 
eastern cases show nearly equal division between cyclonic 
and anticyclonic curvature of the jet stream. 

There appears to be a difference between the western 
and eastern occurrences with regard to the horizontal 
wind shear and the stability. Over 77 percent of the 
western occurrences were associated with a horizontal 
wind shear of less than 25 knots per 100 miles while only 
57 percent of the eastern occurrences were associated 
with horizontal wind shears of less than 25 knots per 100 
miles. Over 46 percent of the western cases had lapse 
rates of potential temperature (A@/Az) greater than 4° C 
per 1,000 feet as compared to only 20 percent of the eastern 
occurrences. 

Figures 3 and 4 relate the horizontal and vertical wind 
shear for the C-A-T occurrences in these two areas. The 
relationship with stability is also indicated through the 
symbol used for the plotted occurrences. Lines indicating 
the criteria of 25 knots per 100 miles for horizontal shear 
and 5 knots per 1,000 feet for vertical shear are drawn 
on the charts. These criteria are close to the 50 knots 
per 150 miles given by Harrison [1] and the 6 knots per 
1,000 feet given by George [2]. 

It is important to note that only 16 percent of the 
eastern C-A-T occurrences (fig. 3), had both the horizontal 
and vertical shear in excess of the above two criteria and 
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Figure 2. Sequence of 300-mb, charts for the April 12-14, 1960 period of pronounced C-A-T activity. Contours (labeled in hundreds of 
feet) are given by solid lines, isotherms (°C.) by dashed lines. Jet axis is shown by arrows. Reports of turbulence are plotted as 
circles (2 light intensity, ® moderate, and @ severe) with height of report in hundreds of feet. 


that 33 percent of the cases had both the horizontal and It is interesting to note that where the vertical wind 
vertical shear below these criteria. However, in 67 per- shear was less than 5 knots per 1,000 feet, 53 percent of 
cent of the cases, either the horizontal or the vertical the eastern C-A-T occurrences had A6/Az equal to or less 
wind shear did exceed these criteria. This suggests that than 1.5° C. per 1,000 feet and only 17 percent had 46/42 
it is not necessary to have both a large horizontal wind equal to or greater than 4.5° C. per 1,000 feet. Where the 


shear and a large vertical wind shear. vertical shear was more than 5 knots per 1,000 feet, only 
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Ficure 3.—Horizontal and vertical shear for the C-A-T occur- 


rences east of 103° W. 


20 percent of the occurrences were associated with the 
less stable situations. This suggests that lower stability 
is necessary in the cases where large vertical shear is 
missing. 

Where the horizontal wind shear was equal to or greater 
than 25 knots per 100 miles, over 58 percent of the eastern 
C-A-T occurrences were associated with low stability 
(A0/Az equal to or less than 1.5° C. per 1,000 feet) and only 
7 percent were associated with high stability (A@/Az equal 
to or greater than 4.5° C. per 1,000 feet). Where the 
horizontal wind shear was less than 25 knots per 100 
miles, the occurrences were associated with a more even 
distribution of the stability values. Apparently, low 
stability is more important in the situations with large 
horizontal wind shear. 

Over 62 percent of the western C-A-T occurrences (fig. 4) 
had both the horizontal wind shear less than 25 knots 
per 100 miles and the vertical shear less than 5 knots per 
1,000 feet. This is in decided contrast to the 33 percent 
of the eastern cases in this category. This would indicate 
that strong vertical 
important over the mountains as over the eastern area. 


horizontal and shear are not as 

Over 33 percent of the western cases were associated 
with low stability (4@/Az equal to or less than 1.5° C. per 
1000 feet) and 37 percent with high stability (A@/Az equal 
to or greater than 4.5° C. per 1000 feet). This is in con- 
trast to the 48 percent with low stability and 13 percent 
with high stability in the eastern C-A-T occurrences. 

This indication of the greater importance of high sta- 
bility over the mountains is in agreement with past work 
on mountain waves [3] which called for a fairly stable 
laver with little evidence of any pronounced horizontal 
wind shear. 

In the April data, there were nine cases which seemed 
to be associated with upper level cold Lows. Six of these 
hine cases were along the west coast. Other C-A-T cases 
seemed to be related to very cold tropopause situations. 
These relationships will require further study. 
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Horizontal and vertical shear for the C-A-T oceur- 


rences west of 103° W 


FIGURE 4. 


5. CONCLUSIONS 


Because of the indefiniteness of many of the original 
pilot reports and the possibilities of incorrect evaluation 
of the meteorological parameters at the scene of the tur- 
bulence, conclusions cannot be drawn with too much 
certainty from the analysis of data in this study. How- 
ever, it is felt that the general indications as stated are 
valid. 

One important fact which must be recognized in such 
studies is that the distance between the present upper-air 
stations is much larger than the scale of the observed 
turbulence. It is unlikely that there will ever be a net- 
work of upper-air stations sufficiently dense and with 
sufficient frequency of observations to insure the determi- 
nation of the exact value of the parameters in the actual 
area and at the time of turbulence occurrence. 

This leads to the suggestion that large-scale patterns 
associated with the development of clear air turbulence 
should be stressed. These patterns can be detected and 
measured on the standard synoptic weather charts. Con- 
siderable attention must then be given to the various 
possible atmospheric processes which will lead to the 
development of these local turbulence situations [4, 5). 

Over the western mountains, apparently both the low 
stability situations which favor mixing and turbulent 
motions and also the more stable wave situations are 
important. East of the mountains, low stability is an 
important factor along with either a strong horizontal or 
a strong vertical wind shear. 

Recent work by Clodman, Morgan, and Ball [6] stressed 
the importance of gravity waves over continental areas. 
They pointed out that the combined gravity wave-vertical 
wind shear mechanism is predominant over land areas. 
The conclusions in this study seem to be in fair agreement 
with their ideas. In spite of the much larger area east 
of 103° W., there were 90 C-A-T cases to the west and 


74 cases to the east of the dividing line. The larger 











YS 


number of occurrences over the mountains supports their 
stress on the influence of the underlying terrain. 

Again, there is the handicap of the lack of knowledge 
about what is taking place in the areas other than those 
reporting turbulence. A more intensive observation pro- 
gram is needed in which reports of all the non-occurrences 
of C-A-T as 
analysis. 

In the meantime, the Weather Bureau is collecting 


well as the occurrences are available for 


weekly a summary of pilot reports covering C-A-T occur- 
rences from all of the Airway Forecast Centers as well as 
various airlines. Further papers covering these studies 


will be forthcoming. 
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THE HURRICANE SEASON OF 1960 


GORDON E. DUNN 


Weather Burea 


1. GENERAL SUMMARY 


The hurricane season of 1960 was subnormal from the 
standpoints of both frequency and intensity. Seven trop- 
ical cyclones developed in the Atlantic and four reached 
hurricane intensity (fig. 1). However, only one hurri- 
cane—Donna—was of major intensity, although Ethel 
may have been so for a few hours. Donna was the most 
destructive hurricane ever to strike Florida and one of 
the most damaging ever to affect the United States. It 
is also believed to have caused hurricane winds over a 
greater proportion of the United States coastline than any 
other known hurricane. ‘Tropical cyclones were well 
distributed throughout the season except in October which 
was remarkably free from even the weak disturbed con- 
ditions normally observed in the Tropics. Tracks of the 
tropical cyclones are shown in figure 1. 

The mean 700-mb. circulation for June [1] was attended 
by above normal heights in the Atlantic over most areas 
north of 30° N. and below normal over most of the 
Tropics and subtropics. This circulation type is usually 
associated with at least normal tropical cyclone activity. 
The June cevelone (unnamed) formed in a manner de- 
scribed by Riehl [2] in which a westerly trough of large 
amplitude extending well into the Tropies fractures and 
the southern portion retrogrades followed by the develop- 
ment of a tropical evelone. 

The maximum positive 700-mb. height anomaly in the 
Atlantic for July was about 150 feet in the middle Atlantic 
between 40° and 50° N. 
trough of considerable amplitude [3] 


with a United States east coast 
Normal tropical 
evelone frequency in the Atlantic area during July is only 
one storm every two vears [4]; thus the two tropical 
evclones represent considerably above normal activity. 
Abby developed in very low latitudes and did not come 
Indeed, 


much of Abby’s track (fig. 1) was along the lowest latitude 


under the influence of the east coast trough. 
of record for July. 

The formation of Brenda was rather similar to that of 
the unnamed June storm. 

Stark [5] noted a major reversal in the large-scale 700- 
mb. mean circulation over North America and adjacent 
ocean areas between the first and last halves of August. 
The belt of positive height anomalies shifted northward 
from 30°-40° N. during the first half to 40°-50° N. dur- 


ing the latter half of the month. Although heights re- 


1 Office, Miami, Fla 


mained near or above normal through most of the Tropics 
and subtropics—a condition not very favorable for tropical 
cyclone formation—the circulation during the latter half 
of the month more nearly conformed to the favorable 
type described by Ballenzweig [6] and Cleo formed on the 
18th. 

In contrast to the subnormal activity in the Atlantic, 
tropical cyclone frequency in the western Pacific in Au- 
gust was unprecedented [5]. Eleven storms were re- 
ported, nine reaching typhoon intensity. <A positive 700- 
mb. height anomaly was located some 1000 miles east of 
Japan with a negative anomaly over Formosa (possibly 
in part due to some six tropical cyclones tracked over or 
near the area) resulting in unusually deep easterlies in the 
formation area. 

Tropical cyclone activity was normal during September 
with two hurricanes, one tropical storm, and several 
tropical disturbances which moved along paths rather 
similar to those of Donna and Florence (fig. 1) but never 
reached storm intensity. The 700-mb. anomaly pattern 
in the Atlantic agreed rather well with the composite 
chart described by Ballenzweig [7j as favorable for tropical 
cyclone development. The mean circulation for Septem- 
ber was predominantly one of high index and the axis of 
mean 700-mb. zonal wind speed maximum was displaced 
a substantial distance north of normal over North Amer- 
ica [8]. 

Tropical cyclone activity in the Atlantic was non- 
The mean circulation for the 
Due to 
blocking, the subtropical ridge and maximum westerly 


existent during October. 
month in the Atlantic area was most unusual [9]. 


winds were displaced to the south of their normal positions 

All but two of the 1960 North Atlantic tropical evclones 
developed in or reached the Gulf of Mexico and all recury- 
atures in the Tropics and subtropics were west of longitude 
70° (fig. 1) indicating an Azores-Bermuda anticyclone of 
considerable constancy and strength June through 
September. 

Damage and fatality statistics are shown in table 1 
Estimates are necessarily approximate. In Donnua’s 
damage total for Florida, citrus damage has not been 
included although 30 percent of the grapefruit and 10 
percent of the oranges were lost. The subsequent increase 
in citrus prices is thought to have compensated for the 
Insurance about 


paid off 
$90,000,000 in property losses in Florida and it is esti- 


crop damage. companies 
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North Atlantic Hurricane Tracking Chart 
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mated about one-third of the property damaged was 


insured weather. 


During the night of June 2: 


INDIVIDUAI 


Tropical: Storm (unnamed), June 22-28.--A_ routine 


REVIEW 


WEATHER BUREAU 


res 


PaBLe | Fatality 
reconnaissance flight into the extreme southwestern Gulf 
of Mexico on June 22 found a large mass of weather 
including heavy thunderstorms and squally winds up to Storn Intensity 
35 kt The lowest observed pressure along the recon- aaveee si 
naissance track was LOOS mb. with no circulation reported, \ Hurricane 
although the Mexican coastal stations did indicate a eer —— 
slight circulation cee... nce 
On the morning of June 23 the 200-mb. level had 
become more favorable for intensification and the barom- 
eter had fallen significantly along the Mexican coast from |... iets 
Tampico to Brownsville indicating a northerly drift of — Flere Storr 
the disturbed condition. A Navy reconnaissance plane 
was dispatched to the area and found maximum winds of — Tet! within United State 


only 15 kt. but sea level pressure of 1006 mb. However, 
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13, 600, 000-4 


and damaae statistics, Vor 
cuctones of 1960 
Dat Damiag Deat 
June 23-26 $3, 600, 000 
July 10 600, 000 ‘ 
July 15 slight 0 
28-29 5, 000, 000 0 
17-20 slight 0 
y-12 3000, 000, 000 13 
4), 500, 000 s 
$0), G00, O00 at] 
13, 000, 000 114 
$00, 000, 000 1ti4 
Sept. 14-1 1, 060, 000 0 
Sept. 22-2¢ slight 0 
$396, 160, 000+ f 


NORTH ATLANTIC ‘TROPICAL CYCLONES 
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NUMBER DATE 
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Poation at 7.00 a m. EST of Date Show 
@ Post * mf 
by ’ 
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eooe® p= 
Underlined pressures in Cleo are estimated as 


it is believed the plane did not fly under the most sever 


24 the tropical storm moved 
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First rain band of tropical storm of June 23, 1960, 
Taken at 


Photo courtesy, Dow Chem- 


showing characteristic curvature Freeport, Tex 
Range markers at 20 statute miles 


ical Co 


inland, south and about 30 miles west of Corpus Christi, 
Tex. Rockport reported sustained winds of 40 m.p.h. 
with gusts to 60 and Padre Island Park 50 m.p.h. with 
gusts to 60. The lowest reported pressure was 1002.4 mb. 
at Alice at 0300 cst on the 24th. 


i tide of 3.5 {t 


Corpus Christi reported 
above mean low water. It appears that 
there was no wall cloud. There was some mild curvature 
on the rain bands seen on airborne radar on the 23d and 
the Dow Chemical! Co. radar at Freeport near noon on the 


=o reported a characteristically curved rain band (fig. 2) 


Apparent cloud centers were reported by radars at 
Victoria and Kelly Air Force Base on several occasions. 
Three 
shrimp boat sank with three lives lost, and another was 
beached. 


attended by heavy rains olf 


fishing piers were wrecked on Copano Bay; one 
The storm moved very slowly on June 24-25 
5 to 15 inches or more from 
Corpus Christi to San Antonio and northeastward with 
considerable flooding Port Lavaca reported 29 76 inches 


of rain for the period June 23-26 Tornadoes were 
reported on the 26th as the dying storm moved north- 
northeastward. 


Unusually heavy rains extended into 


Arkansas and southern Illinois 
Fifteen persons apparently were drowned either in the 


igh seas or subsequent floods Damage, mostly from the 
floods, is estimated at $3,600,000. 
ITurricane Abby, July 9-16 


listurbance which finally grew into hurricane Abby was 


The first indication of the 


received from a ship about 3.5° east of the island of Bar- 
at 0500 Est on July 9. 
reported with east-southeast winds of near 40 kt. 


bac S, Showers weather was 
some 
shower activity had been occurring in the Lesser Antilles, 
and 24-hour pressure changes were small but negative. 


At 9100 gest on the 10th, a report received from the SS 
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Del Oro, located at 13.8° N., 59.7° W., with sea level 
pressure of 1007.6 mb. and wind ENE 45 kt., indicated a 
strong easterly wave or a small vortex. A small center 
passed just to the north of Barbados during the next few 
hours. 

At 0800 est July 10, an advisory was issued on tropical 
storm Abby, based on reports from the Leeward Islands 
and a few ships. The storm was moving toward the 


west-northwest and was forecast to reach hurricane 


intensity during the day. Reconnaissance aircraft were 
dispatched to the area and confirmed the existence of 
hurricane Abby by 1100 Esv. 


mated at 90-100 m.p.h. over a small area near the center 


Highest winds were esti- 


Gale warnings and a hurricane watch were ordered for the 
Virgin Islands and Puerto Rico and for the island of 
Hispaniola as the hurricane moved westward. 

On July 11 and 12, the hurricane continued on a westerly 
course, with doubt concerning its intensity. From recon- 
naissance aircraft and surface ship reports, it appeared to 
be rather poorly organized and much of the time was 
barely discernible on aircraft radar. By the morning of 
the 13th, the hurricane had diminished in intensity, with 
maximum winds estimated at 60 m.p.h. in a few squalls 
near the center in the northern semicircle. By early 
morning of the 14th, the storm had reintensified to hurri- 
cane strength with highest winds of 80 m.p.h. estimated by 
retained this 
strength but remained quite small in size as it skirted along 


reconnaissance aircraft. The hurricane 
the northern cost of Honduras, passing inland in extreme 
southern British Honduras early on the morning of the 
15th. 


moved inland. 


Advisories were discontinued after the evclone 


Considerable rain occurred in Central America from 
Honduras north-northwestward into most of southern 
Mexico and the Gulf of Campeche as the remnants of 
Abby continued west-northwestward over the land area 
Gulf of 


Mexico on the 16th confirmed that the radar center of 


Reconnaissance aircraft in the southwestern 
circulation did not emerge over the Gulf of ( ‘ampeche but 


remained over the rugged terrain of the Isthmus of 
Tehuantepec 
An airborne radat 


picture of Abby (fig. 3), as the 


storm just about attained hurricane intensity off the 


northwestern tip of St. Lucia, shows the characteristic 
figure-nine and the asymmetrical cloud and precipitation 
pattern im the early formative stage at low latitudes 
Only in a few squalls in the main rain band on the north 
and east sides were hurricane winds occurring. The 
radar picture in figure 4 on the morning of July 15 as 
British 
Honduras shows a hurricane with a perfectly round eve 
The eve 


observed by the Swan Island radar as it passed within 


the hurricane was approaching the coast of 


and comptete wall cloud Was apparently not 
approximately 70 n. mi. of the station the day before 

At 500 mb. on July S and 9, the subtropical ndge was 
latitude 27 N. and 
extended eastward into the Atlantic north of the Antilles 


centered over Florida ut about 
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Figure 3 Airborne radar picture of hurricane Abby in formative stages off St. Lucia Note the characteristic figure-9 ech« 
I 
{ ‘ Navy photo 


developing storm Officia 


at about the same latitude. Easterly winds 10 to 20 At 200 mb., beginning on the 9th, moderate to strong 
kt prevailed at Trinidad and northward across the Lee- easterly to southerly flow developed over the easter! 
ward and Windward Islands By 1900 est of the 9th, Caribbean, with the appearance of a vigorous anticvcion 
the winds had become light easterly at Trinidad and off and to the east of the Windward = Islands This 
backed to northeasterly and increased substantially in indicated strong outflow, and favored intensification ol 
the Windward Islands. The winds remained easterly Abby, according to Riehl. On the 12th and 13th, this 
and of moderate speeds in the Antilles until Abby passed circulation had about disappeared although a \ ak 
inland over the ¥ ucatan Peninsula. anticvelone persisted ovel Abby By the morning oO the 

at 500 14th, the antievelone at 200 mb. had again become ore 


Only minor fluctuations were noted in) winds 


mb. at ¢ uracao and San Andres, the direction remaining vigorous and, at the same time, Abby reintensifie Oo 


full hurricane strength 


easterly and the speeds dropping off slightly on the south 
appeared to be fairly well organized whe 


sic ol the storm The subtropical ridge persisted ovel Abby 
the Gulf of Mexico as 


\lexico 


Abby continued westward into first formed and moved into the extreme eastern ¢ rib- 
] 
| 


bean Reconnaissance and surface reports, partici I 
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I RE 4 Hurricane Abby on the morning of July 

cloud Official | 
on the llth, 12th, and 13th, indieated a very small 
circulation, and a minimum of convective activity. This 


is possibly one reason for the loss of intensity. The area 
Was covered with considerable stratified cloudiness, and 
radar coverage Was difficult. 

Peak wind gusts of 55 kt. were observed at St. Lucia 
When Abby passed about over the island, and the lowest 
29.80 Rainfall 


inches mb.) 


pressure was (1009.1 was 
80 inches in 24 hours. Martinique reported gusts 
ot kt 

S 


lives were lost on St. Lucia when the roof of a house 


coll: psed Damage was estimated at $435,000 to prop- 
erty and banana and coconut plantations on St. Lucia. 
\ . . . . 

No «stimate has been received from Martinique, however 


15 as it was approaching the coast of British Honduras 
S Navy Photo 


Note the 


cotnptle te wall 


33 percent of the banana and cane crop was lost, and 


considerable damage resulted to roads and bridges 
Only slight damage occurred elsewhere in the Lesser 
Antilles No loss of life has been reported from Central 
America. Property damage in British Honduras was 


light but damage to crops quite heavy. 

Some remnants of Abby apparently continued across 
southern Mexico and developed into hurricane Celeste 
off the west coast [3]. 

Tropical Storm Brenda, July 28-31 
circulation that can be traced back to a position just off 


A weak cy clonic 


the southwestern Florida coast on July 27 began to deepen 
some 150 miles west of the Tampa Bay area the next day 
By the night of July 28, pressure had dropped to near 1000 
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mb. as the Low began accelerating and moving north- 
eastward onto the Florida coast southwest of Cross City. 
Thereafter the storm continued with gradually accelerating 
speed along the Carolina coasts on July 29, through the 
mid-Atlantic States on the 30th, finally passing through 
the New England States on the $ist and dissipating over 
southeastern Canada 

The storm was not officially named until the 29th when 
reconnaissance aircraft indicated tropical storm structure. 
Karlier aireraft and surface reports indicated rather light 
winds over an area within 50 to 100 miles of the lowest 
A tropical Storm is usually associated with a 
but not until 


pressure 
zone of concentrated winds near the center, 
the Low began accelerating northeastward and had reached 
the coustal area of the Carolinas was this type of pattern 
upparent 

Wind gusts in squalls to 60 m.p.h. were reported from 
many locations along the Atlantic coast and the central 
portion of the Florida Gulf coast. A gust of 65 mi.p.h. 
wus recorded at ¢ upe Cod Canal, however, the highest 
sustained wind at an official Weather Bureau station was 
oS m.p h. at Cape Hatteras. The storm had no oppor- 
tunity to reach hurricane force as the track was mostly 
over land after making landfall on the Florida coast. 

Rainfall was heavy along the entire coast from Florida 
to Maine ranging from over 13 inches around Tampa to 
$f to 6 inches in most other areas. There was considerable 
flooding in the west-central portion of the Florida penin- 
sula. Elsewhere the rainfall was beneficial especially in 
the mid-Atlantic States 

Tides were not excessive due to the rather rapid move- 
ment of the storm and the fact that the center remained 
over or near land While the storm Was developing in the 
eastern Gulf, tides of 3 to 5 feet were reported alone the 
Florida west coast in the Tampa Bay area and south- 
westerly winds in combination with high tides produced 
Waves in excess of 10 feet with considerable erosion along 
the beaches 

Some traffie deaths have been associated with the 
heavy rains, however no deaths can be directly attributed 
to the storm Total property damage was estimated 
near S5,000,000 

Hlurricane Cleo, 
ina broad area of squall) weather some 350 miles north- 
east of Nassau, Bahamas, on August IS. 


considerations, the formation was unique in that a marked 


Auqust 17-20. Hurricane Cleo formed 


From surface 


trough with at least one circulation center was present to 
the north of the area of formation Thus, Cleo’s develop- 
ment did not conform to the normal Increasing easterlies 
and evelonie vorticity in the north portion of the trough 
which usually accompany tropical evclone development. 
Cleo 


parallel to the Atlantic coast at an initial forward speed of 


Hurricane remained small and moved about 


12 kt., later accelerating to 20 to 30 kt. It appears to have 
attained its greatest iIntensits us it approached south- 
eastern New England when winds near the center were 
‘stimated at SO kt 


except during the early and late stages of the storm, 
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there was never a good correlation between reported 
winds and central pressures. Wind reports from aircraft 
were consistently high compared with winds calculated 
from sea level pressures obtained by aircraft penetrations, 
Some compromise has been made and this accounts for 
the estimated sea level pressure noted on the storm track 
(fig. 1). 

There was no known loss of life and very little, if 
any, damage directly attributed to this hurricane. 

Hurricane Donna, Donna, 
the one major hurricane of the season and the most 


A uqust 29-Se pte mber 18. 


destructive ever to strike Florida, was detected by aerial 
reconnaissance on the afternoon of September 2 near 14 


N., 49° W. 


time were 120 kt. 


Maximum observed surface winds at that 
and there was a well-developed eve 
with a central pressure of 973 mb. The existence of a 
tropical storm or hurricane had been indicated earlier in 
the day by surface ships, apparently on the fringes of the 
circulation, which reported shifting winds up to 45 kt., 
and pressures as low as 1004 mb. Although lack of data 
makes it impossible to ascertain the exact time or place 
of increase to hurricane intensity, 1t seems likely that this 
occurred near 40° W. on September 1 and that the early 
stages of Donna were linked to disturbed weather neat 
the African West Coast in late August. 


an active easterly wave through the area Was suggested 


The passage ol 


by unusually heavy rain at Dakar, with which the crash 
of an airliner there on August 29 was associated, and by 
heavy rain in the Cape Verde Islands on the 30th. An 
aircraft also reported indications of a tropical disturbance 


N., 24° W. on that date. Although no additional 


evidence of storm development was received until Sep- 


near 10 


tember 2, evtrapolation of Donna’s track back to thy 
African vives reasonable continuity with these 
disturbed conditions. A study by Dunn [10] shows that 


coust 
tropical evelones of the “Cape Verde” tvpe typically 
behave in this manner, passing through the islands as 
unstable waves and attammime hurricane strength sore 
10° to 15 

On August 29 and 30 a ridge of anomalously-high surface 


to the west. 


pressure extended from near the Azores northeastward, 
as associated with tropical 
Atlantic 
ridge at the 500-mb. level also was unusually strong during 


au pattern long recognized [11] 


evclogenesis in the eastern The subtropical 


this period. However, it extended from northwestern 
Africa to the Azores: a northeastward protrusion toward 


western Europe was not so much in evidence as in some 


other cases of storm formation in this area. At the tim 
Donna Was discovered, a minor tropical depression was 
centered some 500 miles east of Bermuda, moving north- 
ward into a frontal trough. Thereafter, surface pressures 
began to rise to the north and northeast of the hurricane 
and for the next two or three days a ridge of above-nor al 
pressure, surface and aloft, extended westward from just 
south of the Azores, preventing ans marked meridio al 
component of motion. Consequently, Donna continued 
toward the west-northwest on approximately the clima- 


tological track, but at a slightly faster-than-average rite 
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of sbout 17 kt. This course took the hurricane through 
the northern 
September 4 with the eve passing over Barbuda, St. 


Barthelemy, Sint Maarten, Anguila, and about 10 miles 


Leeward 


to the south of Anegada. 
A slight decrease in intensity apparently occurred as 
Islands. The 


maximum sustained wind observed at Sint Maarten was 


the hurricane approached the Leeward 


110 kt. and the lowest barometer reading 952 mb., com- 
pared to earlier reports from reconnaissance of 140-kt. 
winds and dropsonde measurements of 947 mb. Observers 
on the reconnaissance aircraft also had the impression 
that a slight weakening occurred September 5. However, 
no major changes took place during the first few days 
Donna was under observation and the radar presentation 
was characterized as that of an intense, “idealized” 
hurricane, 

Only minor damage was reported at St. Thomas, Virgin 
Islands, with the wind reaching a gust speed of 52 kt. as 
the storm center passed about 35 miles to the northeast 
on September 5. Movement continued toward the west- 
northwest on the 5th and highest sustained winds were 
mly 33 kt. at San Juan, P.R., as the hurricane passed some 
85 miles to the north However, high tides of 4 to 6 feet 
and heavy surf resulted on the northern and eastern 
coasts of Puerto Rico and serious floods developed over 
the northern and eastern portions of the island on the 
morning of September 6. Flood warnings had been issued 
following the detection by the San Juan radar of heavy 
rain moving in on the southern and eastern coastal sections 


1 


on the evening of the dt h However, despite thre warnings, 


107 persons were drowned. The greatest loss of life was 
at Humacao where S4 deaths occurred The people in 


t 


(His area had returned to their homes, many of which 
were built on the river bed, after the hurricane center had 
passed to the north. According to reports, ther failed to 
heed the subsequent flood warnings. The rain totaled 
10 inches or more over a considerable area and measured 
over 15 inches on some of the mountain slopes. 

The track of Donna from September 5 to 7 wes affected 
There 


were temporary indications of recurvature at this time 


by a short-wave trough passing to the north 


but the southern portion of the trough weakened, perhaps 
partially as a result of anticvelogenesis associated with 
the hurricane. Pressures north of the center continued 
to rise as the trough continued eastward and, as a result, 
the hurrieane’s course changed from west-northwest to 
west. In the meantime, Donna had more than regained 
any previous loss of intensity. The central pressure 
given by dropsonde was 940 mb. on the 6th and 944 mb. 
on the 7th. 

Turks Island, in the southeastern Bahamas, escaped 
the ill fury of the hurricane as it passed to the north, 
5O kt. 


fall was extremely heavy and amounted to over 20 inches, 


and highest winds were only 45 However, rain- 


much of which fell in a 12-hour period. Several other 
islands in the southeastern Bahamas were less fortunate. 


Acklins 


lie eye passed over or very near Mayaguana, 


Islands during the evening of 
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No deaths 


were reported but only a few of the more substantial 


Island, Fortune Island, and Ragged Island. 


buildings were left standing in many of the island villages. 
At Mayaguana, which was battered by hurricane-force 
winds for 13 hours, many of the residents took shelter 
in buildings of the missile-tracking base. The maximum 
winds in Donna at this time are not known with cer- 
taintv. At Ragged Island the anemometer failed at 130 
kt. The observer estimated the maximum at Fortune 
Island at 150 kt. after the anemometer and supporting 
tower blew away. Most of the islands in the western 
Bahamas were somewhat north of the storm track and 
no extensive damage resulted although Andros Island 
experienced hurricane winds for several hours. 

On September 9, Donna skirted the northeastern coast 
of Cuba, bringing gales and heavy rains to much of the 
island, then took a west-northwest course, toward the 
Florida Keys. The center crossed over the middle Keys 
just northeast of Marathon between 0200 and 0300 Es? on 
September 10. It is estimated that sustained winds neat 
the center reached 120 kt. and momentary gusts were at 
least 155 kt. The central pressure had continued to drop 
as the hurricane moved across the warm waters of the 
Florida Straits and was approximately 930 mb. when the 
center reached the Keys. Miller {12| has made com 
putations relating the probable minimum pressure = to 
water surface. His 


temperatures of the underlying 


calculations indicate that a temperature of 86° F. would 
be required to produce a central pressure of 930 mb It 
is of interest to note that there are only a lew areas in the 
hurricane belt, mostly in the Gulf of Mexico and neat 
the Florida Kevs, where such high sea-surface tempera 
tures are likely, and that the lowest pressure ever re 
corded in the Western Hemisphere (S92 mb.) occurred i 
the Keys Hurricane of 1935 

Shortly after it passed the longitude of Puerto Rico, 
Donna slowed considerably in forward speed and moved 
at 6 to 10 kt. until it reached the Florida Keys on the 10th 
There were many light and variable winds in the middle 
and upper troposphere and the steering pattern during 
this period Was hever clearly defined At 0700 EsT on 
September 8, when the hurricane was located only 380 
miles southeast of Miami and moving westward at about 
10 kt. the Miami wind at 500 mb. was still blowing from the 
west, It seems evident that there was a contribution from 
the storm itself to its movement 

On September 10 a low pressure system was moving 
across southeastern Canada with a cold front extending 
through New England thence southwestward, becoming 
stationary from Tennessee to Texas No marked surface- 
pressure trough accompanied the front south of about 40 
N. However, at about this time, the ridge of high pres- 
sure which had persisted to the north of the storm began 
to decay. This occurred in response to an upstream 
increase in amplitude in the westerlies and the establish- 
ment of a long-wave trough in the east-central United 
States. Asaresult, Donna began recurvature at this pomt 


and moved northwestward along the southwestern coast 











tadar picture of Donna at its closest position to Miami: 
0730 Est 


igure 5 
illustrating all the features of the classical hurricane. 


sept mber LO, 1960; strobe LOO n. mi. 


of the Florida peninsula. The eye passed over Naples 
Fort Myers as the hurricane turned northward, 
moved inland, and then 
re-enter the Atlantic just north of Daytona Beach about 
0400 Es’, September 11. 
with gusts of 100 to 130 kt. probably occurred on the 


and 
continued northeastward to 
Sustained winds of over 90 kt. 


lower west coast from Everglades City to Naples and 
there were sustained winds of near hurricane force with 
the 
Despite the trajectory over land and a filling 


higher gusts near the center as it moved = across 
peninsula 
of central pressure from 950 mb. at Fort Myers to 970 
mb. on the east coast, the storm was still intense and well 
organized when it moved into the Atlantic again 


The storm surge on the Florida Keys reached as much 


aus 13 feet above normal levels and was generally 8 to 12 
feet some 40 miles to the northeast of the track and 20 


miles to the southwest. Destruction from the combina- 
tion of wind and water in this area ranged from major to 
Tide departures tapered off to the north 


Beach 


On the southwestern Florida coast, the storm surge 


almost complete 
to 3 to 4 feet above normal in the Miami-Palm 
urea 


Wiis locally up to 11 feet above normal, The surge reached 


ur north as Naples, and Fort Myers reported levels 
tto 7 feet On 
Bradenton and the Atlantic coast 


tides were mostly 1 to 3 feet above normal. 


nus 
the Gulf coast north of 
north of Palm Beach, 
Maps showing 


above normal 


the high water marks recovered in the Florida Keys and 
the Florida 
Harris, 


along west coust ol will be presented in a 


forthcoming article, by > i in the Review. 
The effect of the hurricane on the flora and fauna of 


the National 


The largest stand of big mangrove 


Kverglades Park is being studied by spe- 
Clalists in these subjects 
trees in the world is located here and was about 50 per- 
cent or more killed, with a complete kill in some areas 
Almost 
all large mangrove trees which survived the 1935 hurri- 


for reasons which are not vet wholly understood. 
cane were killed by Donna. This is true to a somewhat 
lesser extent of the mahogany trees. It is interesting to 
that mangrove trees growing after the 1935 


1935 


note new 


hurricane were not as large as those standing in 
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indicating a longer period free from major hurricane oc- 
currence prior to 1935 in this area. The great white 
heron, only found in the United States in extreme southern 
Florida and once in danger of extinction, suffered about 
One 


of the two or three largest nesting concentrations of the 


a 35 to 40 percent loss but about 600 birds remain. 


American bald eagle in the United States, exclusive of 
Alaska, is also in the Park. All eagle nests were destroyed 
but four months later some twelve had been rebuilt. 

Figure 5 shows a radar picture at the time Donna was 
closest to Miami. Present are the features of the classical 
hurricane. The round eye about 21 miles in diameter 
and the thick wall cloud some 17 to 20 miles in diameter 
are observable. At the time, the heavy rain band visible 
over the Miami station was producing rain at the rate of 
about one inch per hour. 

Donna began a gradual acceleration as it moved through 
northeastern Florida and was advancing at about 12 kt. 
when it passed into the Atlantic. Rapid re-intensifica- 
tion oceurred over the ocean and when the center was 
about 80 miles southeast of Charleston, 8.C., on the after- 
noon of September 11, the SS \/ae reported winds of 
105 kt. 20- to 30-foot 


tornadic storms were reported in coustal South Carolina 


seas. Several small brief 


and 


with about ten people hospitalized and considerable prop- 


erty damage in the Charleston area. Winds along the 
beaches near Charleston reached about 60 kt. in gusts. 
the North Carolina 


ine received sustained hurricane-force winds. 


However, the coastal section near 

A short-wave trough was moving through the Great 
Lakes region on September 11 and, under the influence of 
an increased southerly flow, Donna accelerated to a speed 
of about 30 kt. 


Carolina coast just northeast of Wilmington during the even- 


toward the northeast, reaching the North 


ing of the 11th, then passing into the Atlantic again near the 
Virginia line about 0500 est of the next morning. During 
its passage over North Carolina, Donna’s eve was unusu- 
ally large with the area of calm or light variable winds 


small- 


ranging from 50 to 80 miles in diameter. Some 
scale irregularities in speed and direction of motion during 
this period Tye be partially attributable to differential 
effect 
suggested as the cause of erratic behavior of some of the 


Minimum 


friction between land and water. Such an was 


1955 hurricanes in this area [13]. pressures 
reported along this section of the track ranged from 958 
to 967 mb. and highest winds were in the 70-90-kt. bracket 
with some estimates of 100-kt. gusts. Tides reached 4 
to 8 feet above normal at various places along the North 
Carolina coast and waves were reported from 15 to 20 feet. 
Several small tornadoes or locally destructive storms 
occurred in the forward portions of the hurricane and at 
least eight persons required hospitalization. Five other 
tropical storms have taken rather similar tracks across 
North Carolina in the past decade and Hardy [14] reported 
that residents in a few areas felt that Donna was the most 
destructive. He points out, however, that this is probably 
true only in limited sections in the northeastern portion of 


the State. 














MivrcH 1961 


When Donna again reached the ocean, it resumed its 
rapid movement with a forward speed of 30 to 35 kt., 
moving northeastward a short distance off the coast and 
crossing Long Island shortly after noon on September 12. 
Winds of about hurricane force, but with gusts locally to 
about 90 kt., brushed the Maryland, Delaware, and 
Residents at Ocean City, 
Ma., described the storm as the worst in the city’s history. 


southern New Jersey coasts, 


Damage to property in other areas along the immediate 
shore was heavy and considerable agricultural losses were 
sulfered inland. 

Sustained winds reached about 90 kt. at several points 
on Long Island and 50 to 60 kt. on western Long Island 
and in New York City. Gusts of 100 kt. or higher were 
reported at Montauk, L.1., and Block Island, R.I., and 
peak gusts reached or exceeded hurricane force east of the 
center through southern New England and northward to 
Winds to the west of the 
track were somewhat less and there were no sustained 


the New Hampshire coast area. 


hurricane-force winds reported on the mainland in New 
England except for isolated cases where local topographic 
effects were responsible. One such exception was a sus- 
tained 80-kt. wind, with gusts of over 120 kt., at Blue Hill 
Observatory at Milton, Mass. The minimum central 
pressure recorded at Brookhaven, L.1., was 961 mb., 
approximately the same as the minimum along the track 
Gradual filling and weakening 


occurred farther north as the center continued rapidly 


through North Carolina. 


northeastward, moving through Maine just west of Cari- 
bou and into Canada late on September 13. Winds of 
hurricane force still persisted in squalls near the center 
until about the time it reached the Canadian border. 
The storm then moved northeastward through Labrador 
and into the Atlantic as a weakening frontal disturbance. 
Complete meteorological data for individual stations in 
the hurricane’s path can be found in Climatological Data, 
National Summary for September 1960. 

There was evidence that Donna was beginning to assume 
extratropical characteristics as early as September 12. 
Many of the maximum winds reported in the Middle 
Atlantic States were from the northwest, indicating that 
cooler, dryer air was beginning to invade the circulation. 
Another feature contributing to peculiar distributions of 
wind and pressure extremes was the unusually large eve of 
During the period the hurricane was moving 
from North Carolina to southern New England this was as 


Donna. 


much as 50 to over 100 miles in diameter, an extreme and 
probably unprecedented size for a hurricane eve 

Tides were well above normal along the entire Middle 
Atlantic coastline and ranged about 6 feet above normal 
near Atlantie City, N.J., and in the New York City-Long 
Islsnd area. Departures of 5 to 10 feet above normal 
occurred along the southern New England coast but 
for‘unately the damage potential was lessened by the 
arr val of the storm surge at the time of normal low tide. 
Even so, considerable damage to coastal properties was 
rted all along the Middle Atlantic and New England 


‘Oust as a result of combined effects of winds and tides and 


re 
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there were extensive damages from flooding farther inland, 
A map 


showing the peak tide and peak surges in New York and 


particularly from the Catskills to Long Island. 


southern New England will be given in a forthcoming 
article, by D. L. Harris, in the Review. 

Donna was unique in that it gave hurricane force winds 
to Florida, the Middle Atlantic States, and New England. 
However, although it was one of the most destructive 
hurricanes of all time, loss of life was remarkably low. 
This can be attributed to three factors timely and 
accurate warnings, effective dissemination by news media 
and other agencies, and the taking of proper precautions 
by the public. The accuracy of the warnings is in large 
part a reflection of the continuous tracking by aircraft 
reconnaissance and land-based radar, which was probably 
the most complete of any hurricane in history. 

Results of the various objective forecast techniques 
ranged from poor to excellent with median 24-hour errors 
for the individual methods ranging from 68 to 132 n. mi. 
The results of the various systems cannot be compared 
since the number of forecasts and the times of the fore- 
casts with each were not the same, being dependent on the 
available data. Each method, including the one with the 
lowest median error, showed serious errors during at least 
one of the three critical periods the abortive recurvature 
off Puerto Rico, the recurvature over Florida, or the 
acceleration off the Atlantic coast. 

Hurricane Ethel, September 14-17. 
developed rapidly in the central Gulf of Mexico early on 
September 14. Its position and intensity were established 
by the 0930 est report from MAMOS (Marine Automatic 
Meteorological Observing Station) in the central Gulf of 


Hurricane Ethel 


Mexico. The hurricane moved northward and continued 
to intensify rapidly during the day with a central pressure 
of 972 mb. and winds of 140 kt. reported by reconnats- 
sance aircraft that afternoon. During the night of Septem- 
ber 14-15, cool dry air entered the circulation and the 
hurricane’s intensity diminished quickly. The hurricane 
center reached the coast near Biloxi, Miss., with the lowest 
pressure 981.4 mb. during the afternoon of September 15 
at Keesler Air Force Base. 
moved northward through eastern Mississippi that night 


It continued to weaken as it 


The remnants of the storm were located in central 
Tennessee on the morning of September 17. 

The highest sustained wind reported by a land station 
was 78 kt. with gusts to 90 at Venice, La., at 0415 csr, 
September 15. Burrwood, La., reported winds of 45 kt 
with gusts to 60. The highest tide reported was 7 feet 
above mean sea level on Quarantine Bay on the east side 
of the Mississippi River about 0400 csr on the 15th 
Tides of 2 to 5 feet above mean sea level occurred from the 
mouth of the Mississippi River eastward to St. Marks, Fla 
Local rains of 5 to 6 inches or more were reported through 
southeastern Mississippi, southwestern Alabama, and ex- 
treme northwestern Florida but no significant flooding 
occurred. 

The total storm damage probably did not much exceed 


one million dollars and there Was ho loss of life. 
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Tropical Storm Florence, September 17-26.—An extensive 
shower area was noted well to the northeast of the Leeward 
Islands September 16 although there were no indications 
of a definite circulation. By the morning of the 17th, 
pressures through the eastern Antilles had fallen 3 to 5 mb. 
with light south and southwest winds indicating the possi- 
bility that a circulation had developed. On the evening of 
the 17th, reports from shipping to the north of Puerto Rico 
placed a closed circulation near 21° N., 66° W. 


winds up to 35 m.p.h. The Low continued westward 


with 
ubout 10 m.p.h. and reconnaissance aircraft located a 
broad, ill-defined center with maximum winds around 
10 m.p.h. on the morning of the Sth near 21° N., 69° W. 

Tropical Storm Florence moved on a west to west- 
northwest track near 12 m.p.h., gradually weakening until 
the 20th when reconnaissance aircraft found only a weak 
Low south of Andros Island in the Bahamas with no 
significant weather or strong winds. The remains of 
Florence moved into the western end of Cuba and became 
nearly stationary until the evening of the 22d when con- 
ditions became more favorable for redevelopment. The 
Low began moving northeastward and was located just 
off the southwestern Florida coast by the morning of the 
23d with winds up to 30 in.p.h. and widespread rain over 
southeastern Florida. It then became blocked by a large 
high pressure system along the mid-Atlantic coast after 
reaching the vicinity of Lake Okeechobee the evening of 
the 23d and changed to a west-northwest track, drifting 
into the eastern Gulf of Mexico near Tampa early on the 
25th 


moved into the Pensacola area on the morning of the 26th 


The Low continued quite weak over the Gulf and 
with winds less than 25 m.p.h. but with a rather large rain 
area that covered the southern portions of Alabama and 
Georgia and northwestern Florida. 

Florence was never a well-defined tropical storm and 
aximum winds were just barely of tropical storm in- 


(for only a short period) although gusts to 52 


tensity 
in.p.h. were reported in the Vero Beach area in a squall 
when the Low was nearest that station. The only signifi- 
cant damage in this storm was from local flooding in 


Florida. Rainfall totals of 3 to 6 inches or more were 
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reported during passage of the Low on ground that was 
already saturated from the previous heavy rains of Donna. 
Monthly totals for September were well over 20 inches at 
many spots in southeastern Florida with a few totals in 
excess of 30 inches. No damage of consequence was 
reported over the eastern Gulf States with rainfall prob- 
ably beneficial to crops. No known fatalities or injuries 
have been reported. 
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THE WEATHER AND CIRCULATION OF DECEMBER 1960' 
An Unusually Cold Month in the United States 


ROBERT H. GELHARD 


Extended Forecast Section, U.S. Wea 


1. HIGHLIGHTS 


Average temperatures throughout most of the contigu- 
ous United States during December 1960 underwent a 
reversal from the regime which had 


marked warm 


dominated the fall season [1, 2]. Subzero temperatures 
were frequent in the northern half of the country east of 
the Rocky Mountains, and freezing temperatures threat- 
ened winter vegetable and citrus crops from southern 
California eastward through the Rio Grande Valley to 
central Florida. 

Snow cover, which had been restricted to the eastern 
slopes of the Rockies and extreme northern Plains, spread 
throughout the northern half of the country early in 
December and for the most part remained throughout 
the month. Inthe southern Plains, precipitation in excess 
of normal resulted in some flooding around mid-month, 
while drought conditions in the Great Basin which had 
been alleviated somewhat in November resumed during 
December. 


2. THE GENERAL CIRCULATION 


The general circulation at 700 mb. for December 1960 
fig. 1) represented a distinct change from the predomi- 
nantly zonal or high-index state of the preceding month 
I], particularly in the western portion of the Northern 
Hemisphere from western Europe to the central Pacific 
Ocean. 


height pattern for the preceding month shows that the 


A comparison of figure 1 with the mean 700-mb. 


diffuse planetary wave system of November became re- 
solved into a series of well-defined troughs and ridges 
located in climatologically preferred areas [3]. 

The field of anomalous 700-mb. height changes (changes 
after the normal change has been removed) from Novem- 
ber to December (fig. 2) helps delineate the areas of great- 
est amplification. The greatest change, +680 feet, took 
place in the central Atlantic with the development of the 
Azores High at a higher than normal latitude. This was 
associated with a northward shift of the mean westerlies 
and subsequent deepening of the downstream trough in 
the eastern Mediterranean where anomalous heights fell 
by 310 feet (fig. 2). 


show that amplification of the wave pattern of the gen- 


The 5-day mean charts (not shown) 





A les describing the weather of January, February, and March 1961 will appear 
April, May, and June 1961 issues, respectively, of the Monthly Weather Rericu 


eral circulation apparently began with the rather sudden 
development of the Azores High late in November. Sub- 
sequently this amplification seems to have spread up- 
stream, first to western North America and then to 
northeastern Asia. 

Relatively little change occurred in the eastern part of 


Nearly all 


of Asia was dominated by a vast Siberian High at sea 


the hemisphere except in northeastern Asia. 


level, extending from the Black Sea to Japan (fig. 3). 
Generally stormy conditions existed in Europe and the 
Mediterranean area, but the modifying effects of maritime 


air did not penetrate bevond European Russia. 
3. THE INDEX CYCLE 


The change in the westerly circulation from November 
to December was a classic example of an index cycle as 


described by Namias [4, 5]. Figure 4 shows the march 


»- --0O 


of indices of temperate (35°-55° N.) and subtropical 


(20°-35° N.) westerlies at 700 mb., as defined by over- 
lapping 5-day mean maps computed three times per week. 
The index referred to here is for the western half of the 
Northern Hemisphere from 0° to 180° W. 
westerlies reached a peak near the end of the third week 


The temperate 


in November and then began a four-week decline before 
The dotted 
line, giving the daily components of the mean index, is 


making a rapid recovery late in December. 


included to show the very sharp drop in the index which 
is damped somewhat by the averaging process and coin- 
cided with the first of several severe storms in the South- 
west. The decline in the temperate latitude index was 
accompanied by southward displacement of the westerlies 
except in the North Atlantic, so that the subtropical 
index rose as rapidly as the zonal index fell, reaching a 
peak December 21, a day after the minimum in the 
temperate index. 

Although there seems to be a strong preference for 
index eveles to occur in late winter, they are not uniquely 
fixed with respect to month. One necessary condition for 
the onset of an index evcle seems to be an abundant 
supply of cold air in polar regions. Figure 5A shows that 
the 1000—-700-mb. thickness for November averaged 130 
feet below normal in northwestern North America, thereby 
satisfying this condition. However, in order to maintain 
the atmospheric heat balance, the containment of polar 
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FIGURE | 


Mean 700-mb. height contours (solid) and departures from normal (dotted), both in tens of feet, for December 1960 


Stronger 


than normal northwesterly flow between the intense ridge in western North America and the deep trough along the east coast produced 


very cold weather in most of the contiguous United States 


air brought about by fast westerlies must end at some 
point. This point was reached in December, resulting in 
completion of the index evele and redistribution of cold 
and warm air masses, as strikingly illustrated by the 1000 
700-mb. thickness anomaly patterns for December and 
November, figures 5A and 5B, showing a complete reversal! 
The effects of this 
exchange on the surface weather will be discussed in more 
detail below 


in sign of the thickness anomalies. 


4. TEMPERATURE 


The decline of the zonal index and amplification of the 
circulation pattern through the first three weeks of De- 
cember had a marked influence on the weather of the 
United States. Persistent northerly anomalous flow east 
of the Continental Divide (fig. 1), closely paralleling th 
sea level isobars (fig. 3), brought colder weather to 10s! 


of the contiguous United States. Of a hundred stations 
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NOV. TO DEC. 1960 











Figure 2.—Difference between monthly 700-mb. height 
anomalies for November and December 1960 (December minus 


November) in tens of feet. 


mean 


Anomalous rises in western North 
America and falls in the east were associated with amplification 
of the long-wave pattern. 
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Figure 3.—Average sea level pressure (solid lines) and its de- 


parture from normal (dotted lines) in millibars for December 
1960. High pressure was accompanied by cold weather in the 
contiguous United States. 


distributed uniformly across the nation, 78 percent de- 


freased in average temperature by at least one temperature 


class (out of five classes) from November to December 
fig. 6A). Of the remaining 22 percent, only 3 percent 
increased. Even in the area of no change or increase, 
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NOVEMBER DECEMBER 

Figure 4 700-mb 
meters per second for the western zone of the Northern Hemi- 
sphere from 0° to 180° W. 
35° and 55° N., and subtropical index between 20° and 35 
N Values are for 5-day 
The dashed line represents a daily index, and the dotted lines are 
the normals based on the data of Weather Bureau Technical 
Paper Vo | 


Time variation of 5-day mean indices in 


Temperate index applies between 


periods ending on dates indicated 


only in the extreme northern Plains were temperatures 
above normal for December. 
The Great 


under the influence of large stagnant anticyclones. In 


Basin area and California remained cool 
spite of abundant sunshine, record daily minimum tem- 
peratures in the cool dry air mass kept mean temperatures 
below normal. An interesting effect of these persistent 
Basin Highs was reported from Boise, Idaho, where the 
average wind speed was only 4.6 m.p.h., the lightest in 
21 vears. 

The most unseasonably cold part of the country was 
the East and South as far west as Arizona (fig. 7). Figure 
8A shows the trajectory of migratory anticyclones carry- 
ing cold air masses from Canada just west of the Great 
Lakes and through the Ohio Valley. Throughout the area 
from the Mississippi Valley eastward some record daily 
minimum temperatures were established along with record 
9° F 


low temperatures for the month, such as at three 














































Mean departures from normal of L000—700-mb. thick- 
November B 1960, 
Note the reversal in phase as the index 
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FIGURE 6 \) The number of classes the anomaly of mon 

mean temperature changed from November to December 1 
and (B) the number of classes the total precipitation category 
changed from November to December 1960 Areas of positive 
change greater than two classes are hatehed, and areas of negative 

12 change greater than two classes are stippled. 
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at Evansville, Ind. 
For most stations, average December temperatures were 
10 


In Florida, freezing temperatures threatened 


stations in Kentucky and —7.5° F. 





the lowest in from to 25 vears and close to record 


mtensity 
citrus and vegetable crops as far south as Tampa on five 
Richmond, Va. reported December 1960 as 


Occuslons 


the only month on record with every daily minimum 














temperature below 32° F. 
In the lower Mississippi Valley and Southern Plains _ fF. 
anticvelonic conditions prevailed at sea level (fig. 3), ‘ a |] : . 
emanating from Great However, : , . 





chiefly Basin Highs. 
there was also some contribution from polar air masses, as oe 
for De 


associated 


by 
ber 


Temperature departure from normal I 
[6] Widespread cold 
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evidenced by subfreezing temperatures in the Rio Grande 1960 (from ie 


Vallev and resultant crop damage circulation shown in figures 
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Below normal temperatures in Arizona and New Mexico 
were maintained by persistent cloudiness due to cyclonic 
conditions aloft (fig. 1) which were quite pronounced dur- 
ing the first three weeks of the period. In addition, heavy 
snows early in the period probably contributed to lower 
than normal temperatures. Late in the month a warm- 
ing trend began in the Southwest, coincident with the 
rising zonal index. 

While new records for daily minimum temperature were 
common throughout most of the country, warm air in 
advance of a severe storm in the central Great Plains pro- 
duced record high maximum temperatures in the North- 
east near the end of the first week in December. 

The amplification of the trough-ridge system in the 
eastern Pacific and western North America, which pro- 
vided the mechanism for bringing cold air into the con- 
tiguous United States, was also responsible for the warmest 
December on record for most of Alaska. A strong anom- 
alous southerly component of the average flow, shown by 
the height anomaly pattern in figure 1 (coincident with 
the sea level isobars, fig. 3), produced an influx of warm 
maritime air throughout most of the month. The extent 
of the exchange of air masses is well illustrated by the 
change in the 1000-700-mb. thickness charts for November 


and December (fig. 5). 
5. PRECIPITATION 


In general, excessive precipitation in the contiguous 
United States was restricted to the central third of the 
country and to the North Atlantic coast (fig. 9). An un- 
usual feature of December’s precipitation regime was the 
large amount of snowfall occurring so early in the season. 
In the Southwest and the Northeast many new monthly 
snowfall records were established, such as 10 inches at El 
Paso, Tex., and Winslow, Ariz., 17 inches at Dayton, 
Ohio, and 21 Worcester, Mass. Moderate 
amounts were reported elsewhere except for the Great 


inches at 


Lakes and most of the Great Basin where representative 
stations such as Lansing, Mich., and Salt Lake City, Utah, 
reported their driest Decembers on record. In contrast, 
Waco and Dallas, Tex., reported the wettest December 
on record. 

It is interesting to compare the precipitation patterns 
for November [1] and December in terms of high and low 
index circulations. _November’s pattern was wet in the 
west and relatively dry east of the Rocky Mountains—a 
pattern often associated with fast westerly flow. How- 
ever, in December amplification of the ridge in western 
North America and deepening of the upstream trough pro- 
duced a strongly meridional steering current aloft, and the 
westerlies were weaker than normal throughout the West 
fig. 10). 
enough to block effectively the passage of migratory 
evclones from the west (fig. 8B), so that the Pacific North- 


west during this low-index regime received from 25 to 50 


The intensity of the western ridge was great 


percent less precipitation than normal, and the Great 


Basin received less than half of normal. On the other 
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Ficure 8.—-Number of (A) antiecyclone passages and (B) cyclone 
passages (within equal-area Guadrilaterals of 66,000 nn. mi 

during December 1960. Primary tracks of migratory sea level 
systems are indicated by solid arrows Areas of zero frequency 


are stippled 




















Ficure 9 Percentage of normal precipitation for December 1060 
from [6] 


the Mississippi River and the Continental Divide but deficient 


Amounts were generally in excess of normal between 


elsewhere 
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Mean 700-mb. isotachs in meters per second for 
December 1960. Solid arrows indicate primary axes and dashed 
arrows secondary axes of maximum wind speed 
from normal of mean wind speeds for December 1960 (in meters 
weather in the Great weather west of 


Continental Divide 


hand, the weakness of the westerlies diminished foehn 
drying and favored upslope precipitation in the Plaing 
region between the Mississippi River and the Continental 
Divide (fig. 10B). 

Another common feature of low-index circulations is 
the presence of areas of mid-tropospheric confluent flow 
at fairly low latitudes such as that found in the Texas 
area (fig. 1). While this was not an unusually strong 
confluent zone, there was sufficient contrast between cold 
and warm air masses to lead to three major storms in 
the southern Plains and up to 400 percent of normal 
precipitation with some flooding in eastern Texas and 
southern Arkansas. These same storms produced heavy 
snows as they moved through the Ohio Valley, roughly 
along the track shown in figure 8B. They were responsible 
for record December snowfalls from northern Virginia 
through New England as they deepened in the area east 
of the mean trough near the Middle Atlantic coast (fig. 1), 

Figure 6B illustrates well the general change from the 
precipitation regime of November to that of December; 
that is, from wet to dry in the West, a general increase in 
the Plains States, and little change elsewhere. While it 
has been mentioned that November had a typical high- 
index precipitation pattern, December's pattern was not 
necessarily typical of low index since distribution of pre- 


cipitation is dependent on the precise geographical loca- 


tion of the quasi-stationary trough and ridge systems, 
which can exhibit considerable variation during p 
of low index. 
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